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SUMMARY 


This report covers work performed during Phase I of a four phase 
aggressive analytic/development program whose goal is to reduce 
helicopter passenger noise exposure to the levels experienced in 
fixed-wing coiruflercial jet aircraft. The following is a brief 
summary of the Phase I effort. 

A review of information relevant to rotorcraft interior noise was 
performed with emphasis on the dominant acoustic and vibration 
noise sources and the transmission paths responsible for noise in 
the passenger compartment (cabin). A rank ordering of these noise 
sources established the main rotor gearbox as the primary contri- 

noise levels. In addition, current noise control 
methods were found inadequate for development of general noise 
control concepts . 


The major portion of this Phase I effort involved the development 
ot a detailed analytic modeling approach with capabilities for 
evaluating a variety of noise control concepts (exclusive of 
source modification). This required a model containing details of 
transfer from various sources, via a multitude of 
paths, to the cabin acoustic space. Since the number of struc- 
tural and acoustic resonances in the frequency range of interest 
s^^'tistical approach (Statistical Energy Analysis - or 
SEA as It IS commonly called) was applied. Included in the model 
are many airframe parameters, such as coupling and damping loss 
lactors and source attachment point impedances. 

A comprehensive measurement program was devised to validate the 
entire model during the Phase II effort. Ground test measurements 
will be used to evaluate the airframe parameters and energy path 
transfer functions. Flight test measurements will determine 
source levels for input to the model. 
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INTRODUCTION 

The control of interior noise is a continuing problem in aircraft 
particularly rotary-wing aircraft. The interior noise levels i A 
otorcraft are higher than in fixed wing aircraft because the 

^^omponents are in close proximity to the cabin 
nd are rigidly connected to the airframe. The higher noise 
Speech interference and increased voJLl output 

fnd if for conununication. The result is a decrease in passenger 
and crew comfort. 

I!!d sources in rotorcraft are the main rotor gearbox 

systems; while the engines, turbulent boundary 
fEcm' rotor, and Environmental Control Units 
/Ik sources. Figure 1 shows schematically these 

sources and the paths to the cabin acoustic space. 

gearbox, vibratory excitations are produced as part of 
process and these excitations result in vibration 
housing and airframe which ultimately radiate 
ai the Cabin. Hydraulic system noise exists in some 

such as the s-76, and is a function of pump selection and 
attachment of the hydraulic lines to the structure. 

annarent aerodynamic origin and becomes 

apparent at high flight speeds. The boundary layer noise is 

cockpit windshield and 

sufficient transmission loss to attenuate the boundary noise 
Engine noise depends on the choice of the power plant and Its 
location relative to the cabin. Casing radiated noise, inlet, and 

ECU^noiQ^^s®® coupled to the cabin via airborne paths. 

noise IS generated by the air-cycle machine or blowers and is 
generally conducted to the cabin via ducts. Noise Ly also be 

into^the^cabin^^^ grills at the point where the air is exhausted 

fhfi- bhe cabin, which is typical of 
that achieved in fixed wing commercial aircraft, each of the above 

subjected to noise control. Literature revtefhls 
control measures for rotorcraft have been 

havA^h^^.n noise prediction techniques so far 

limited to empirical relationships that are valid only 

dAVAifIloH which they were derived. The procedures 

^ predicting noise levels inside the cabin cannot 

handle certain types of noise control measures such as vibration 
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structural damping. m most iLtancel airframe, and 

DeFelice [l]) the complex orobliL^?,. example: Levine and 

paths is ac.t'.nowledged but thiSe structureborne noise 

with a uniform intensity distr^l^^ are then quickly replaced 

surfaces. Thus limitations ier^ni I dominant radiating 
involving structural and source/nai-h^^m^^^f trade-off studies 
technology now allow for modelinq of Advances in 

these complex helicopter airframe ^structure! vibration in 

thLe'^advance^LcLun^^ application of 

airframe and tL iTte^lr spa°ce'' 

-loise levels, refinement of existina - ^^^leve reduced cabin 

as improved predictive me^hr>Ho noise control methods as well 

approach to rld^cing^IL To^se Additionally, ?ie 

Since past experience has indicated thai- a systems basis 

one noise source may enhance the nois^^ ^ alleviate 

It IS necessary that the methods hf another source. Also, 

basis to insure that developed noi< 5 #a ^ ®*^alytical 

beneficial to a wide rana^ control techniques will be 

incorporated in the early Itaae^ and can be 

trade-off studies. ^ ^ aircraft design for economic 

noise" redict\on p?ob?L wit^^res^'p'e^c^^^^^^ 
etc., since other research ef fortes ^ar^^^ a 

Rather the intent is to focus on -ht m this area, 

paths exterior to the gearbox or otr^o- '^^^^^^^^y/^coustic energy 
the cabin, to modify the paths of through and into 

develop vihratron isolation systems to%e"Sce IJf ct\n“noise 

strate minimum wei^t ^ectoigues^o¥^'?BH^^ a"< 3 demon- 

noise levels of a rotorcraTt 1,Tth a aoa? k*"® interior cabin 
tive equivalent noise levels of curr^nV «? subjec- 
The emphasis of the present task is to * 

Sion and other noise control ® develop vibration suppres- 

of the gearbox/aiJ??aSe%abL spa^^^ the patL 

of the physical and operational within the framework 

craft and its flight en^f^o^ent a ® ^<^tor- 

concepts shall be developed ^d validatid k ® control 

order to understand the Lsic^obtli^t^ ^ experimental tests in 
apply the techniques to%uturl imirn^'s^.^ mechanism involved and to 
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Tail 8 Main 
Rotors 


Figure 1. Internal Noise Sources and F^aths 



DESCRIPTION OF HELICOPTER NOISE 


Commercial Helicopters in General 


current co^!r^ial h^'^Uc^ters'' are^ 

nyarau.ics, and engine noise sources [3, 4 5 6 7 fti tk?' 

major helicopter manufactuiers (both na?f Infl '’a^rinterait^SL^* 
govIrim“n\s®''«ar?d/^^^^^ amongst manufacturers, operators, and 

considIration^ "'°ihis''®car^°”l othe^J^ade-off 

cabin acoustic properties. caoin, ana the 


Noise Sources for Two Specific Helicopters 
and r^°Ltra“iror,?^ '"I"’ «« ^^own in Figures ; 

SF-rsSS 

Stiff relative to th^ s^in which is quite 
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ge^ailed case study o f Agusta A-109 noise sources. - The Agusta 
A-109 IS designed to provide high performance rotary-wing aircraft 
tor the business and commercial aviation market. The manufacturer 
recognized that some reduction in cabin noise was necessary to 
provide maximum customer acceptance. Therefore, the company 
embarked on a noise reduction R&D program in 1976. A flow chart 
fZZ program is shown in Figure 4. Figures and information 

tor this section are based on data published in Reference [8]. 
The goal of the program was to reduce the cabin noise to a Speech 

5 a 2 > 84 dB for a -Quiet- interior 

and 76 dB for a "Silent'^ Ihtfetiir. An A-weighted sound pressure 
level of 83 dBA was set as a target for the silent interior. For 
ypical cabin noise spectra/ achievement of the 83 dBA taxaet also 
achieves the 76 dB PSIL goal. 


The goals of the noise reduction program were substantially met 
without structural modifications. The resulting treatment is 
commonly referred to as a "carpet hanging" treatment and includes: 
(a) covering all non-window surfaces above the seat cushions with 
plastic trim or upholstery, (b) use of "double wall" construction 
P^bsis, (c) use of vibration isolation to prevent vibra- 
tion transmission to the trim panels and air conditioning/heating 
ducts and additionally on the "Silent" interior, (d) use of a 
clear plastic divider between the pilot and passenger compart- 
merits / and (e) use of double windows in the passenger compart- 
ments. In addition to these basic carpet hanging treatments, it 
was necessary to improve door gasketing and use heavier door 
hinges, and to use air conditioning/heating unit silencers. 


The weight of the overall "Silent" treatment is approximately 115 
^ one-third octave band spectrum for the cabin noise in 
the Silent" interior is shown in Figure 5 and compared to levels 
on the stripped interior. 
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Figure from Reference 8 
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reauction program was to identify ^the^ 

tion/sound transmission paths ^La4 th?/® and vibra- 

cabin noise were taken for a stripped" Lt^for fs'j of 


The correspondeLe of li^es in thi ioi e Figure 6. 
and harmonic frequenciel of i-h^ h spectrum with fundamental 
serves to identif^tSet^ comLnenf rotating machinery 
primary source is^he main gJ^rbox TaiT^ro^^ sources. The 
cant in the frequency range^from lOO to Jnn ^ signifi- 
component of the A-weighted noise^°°Th ® "*ajor 
butes to the spectrum belL lo^ Hz contri- 
source of A-weighted noise Althouah tht i ^ ignored as a 
spectrum make the major coAtribStioo^i-^^fv,^^^®® narrow band 
level, the broadband Lu^^el are Si - ^ A-weighted cabin noise 
are associated with the fluctuating o "^^^se sources 
layer , broadband components of roto^ boundary 
in the turbines, it is difficult t^n ciJf ®^' combustion noise 
butions of these broa^and sourcL "^^e relative contri- 


in Figure 7. Note that peaks in the one^th^^rt^^ speeds are shown 
800, 1600, and 4000 Hz corresnond to octave spectrum at 
spectrum and dominate the overall levels narrow band 
band spectra for the different «nooHo Analysis of the narrow 
ponents decrease with increasina line corn- 
components increase with incre^Iina flLff®®'^ broadband 
acting effects are not unexoect/^^ ^^®®® counter- 
fication of the relative st^na^i of tJ!f k ® ^^^ntitative identi- 
sible. ^ ^ the broadband sources impos- 


dominant sourc^e°^j^n ^this^ catec^^ identified for the A-109. The 

cabin heating. category was the compressor bleed-air for 





SOUND PRESSURE LEVEL IN dB re 2x 10 N/m 



f^igure 6 Narrow Bandwidth (6.2 Hz) Spectrum of A-109A 
Cabin Noise Stripped Interior - Speed 72 ni/s 


Figure from Reference 8 
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Figure 7 


)2S 250 500 1000 2000 4000 gOOO 16.000 32.000 

ONE-THIRD OCTAVE BAND CENTER FREQUENCIES IN Hz. (cps) 


SPL vs Speed - Center of A-109 Cabin Stripped Interior 


Figure from Reference 8 






Noise transmission paths . - The correspondence of the lines in the 
narrow band spectrum with harmonics of the rotational frequencies 
serve to identify the sources of noise, but do not indicate the 
path by which the noise is transmitted to the cabin. The path 
identification for the A-109 was carried out in two steps. Since 
the objective of the noise reduction program was to reduce cabin 
noise by treatment of the radiating surfaces in the cabin, the 
first step was to obtain average vibration levels for each surface 
in order to identify their relative contribution to the cabin 
noise . 

The time-average acoustic power radiated by a surface can be 
written as: 

'"^rad = <Po<=o^ A ‘^rad 

where W is the radiated power, p c is the characteristic 
impedanc^ of the cabin acoustic space,® A is the surface area of 
the radiating surface, <v^> is the mean-square vibration averaged 
over the surface, and u ^ is the radiation efficiency. If a , 
is the same for all sinntaces, then the area-weighted vibratfSn 
level Co be used to rank order the radiating surfaces. 

Measured results are shown in Figures 8a and 8b for the stripped 
interior. The major sources are the aft bulkhead and the overhead 
panel. The windows are not a major source of noise in the strip- 
ped interior. However, they can become a major source when the 
other radiating surfaces are treated to reduce their noise radia- 
tion. 
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Relative Panel Contribution to A-109 Cabin SPL - 72 ni/s, 
Strioped Interior (Area Weighted Panel Vibration) 


Figure from Reference 8 



The neesured contribution to the cebin noise by the various 
surfaces can be used to calculate the noise reduction needed to 
meet a given overall objective. This is done in Figures 9a and 
9b. The procedure used was to assign a noise reduction target for 
each surface such that the contribution by each surface to the 
ov®^^ll noise is egual. Note that the goal of 83 dBA reguires 
that the noise radiated by the windows be reduced. 

The second step in the path identification was to determine the 
contribution of various source/transmission path combinations to 
the cabin noise. A transfer function technique was used. In this 
technique measurements were carried out on the ground using 
artificial sources - loudspeakers and shakers. The artificial 
source is used to excite the structure at a point corresponding to 
an input location for one of the sources. The cabin noise is then 
measured and the ratio of rms noise to rms input vibration is 
obtained as a measured transfer function. For example, to measure 
the contribution by airborne gearbox noise, a loudspeaker is used 
to excite the acoustic space surrounding the gearbox and the 
acoustic level in the space is measured as the input to the 
transfer function. The resulting noise at several points in the 
cabin is measured as the output and a transfer function is deter- 
mined as a function of frequency. To measure the contribution by 
structureborne noise from the gearbox, a shaker is used to excite 
the frame at a point where the gearbox attaches. The vibration 
leve_ is measured at the attachment point as an input to the 
transfer function and the cabin acoustic level is measured as an 
output. The process can be repeated to determine transfer func- 
tions for each gearbox attachment point. 

After the transfer functions have been determined in a ground 
test, measurements are made in flight to determine source levels. 
The rms source level is multipled by the approximate transfer 
function to determine the contribution of the particular source/ 
path combination being considered. The contribution from all 
sources and paths is suinmed incoherently (addition of mean-square 
sound pressures) to obtain the overall sound pressure level. This 
summed level can then be compared with the measured overall level 
to determine the adequacy of the transfer function analysis. If 
important source/path combinations are left out, the measured 
level will be significantly higher than the predicted level. 
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(83 dBA Based on .5, 1, 2 kHz) (Noise Reductions are 
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Figure 9B. Reduction of Sound Radiated from Major Surfaces Required 
to Achieve Target Noise Level for "Silent" Interior (83 
dBA Based on ,5, 1, 2 kHz) (Noise Reductions are 
Referenced to Stripped Interior) 


Figure from Reference 8 
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DIAGNOSTIC TECHNIQUES FOR SOURCE AND PATH IDENTIFICATION 
Review of Existing Case Studies 

The techniques for source and path identification used in the 
existing case studies have been predominantly experimental in 
nature. Measurements are performed in flight for different 
operating conditions . Acoustic and vibration sensors are placed 
at appropriate locations near the sources, along the transmission 
paths, on the interior surfaces of the cabin, and at passenger 
locations. The selection of sensor locations is based on past 
experience and understanding of the cabin noise environment. 

Source identification is often crudely based on the highly tonal 
character of the cabin noise spectrum. The gearbox is an im- 
portant source of distinct tones where the tonal frequencies are 
straightforwardly related to gear tooth mesh rates. Other mech- 
anical occurrences producing distinct tones include turbine blade 
passage, shaft rotation, rotor blade rate, hydraulic pump piston 
rate, etc. In identifying sources the observed tones in the cabin 
noise spectrum are matched with the known mechanical processes 
that produce tones. It is not possible to do the same for broad- 
band sources, neither is it possible to identify the transmission 
path(s) for such sources once they have been identified. 

When the source contributions are not distinctly identifiable in 
the cabin noise spectrum, then the identification is often made 
based on a relatively simple modeling of the system. Source 
levels are either measured directly or estimated analytically and 
combined with a transfer function for the transmission path in 
estimating the contribution to the cabin noise level. 

Measured acoustic source levels are generally combined with a 
transmission loss model for the intervening skin panel, window, or 
bulkhead separating the source from the cabin interior. Differing 
levels of sophistication are used in characterizing the source 
based on its spatial characteristics [9]. Often an attempt is 
made to characterize the directivity of the main or tail rotor 
noise in quantifying the exterior levels over the cad)in skin 
panels. The spatial correlation of the exterior pressure field 
exciting the skin panels is also important in estimating the 
acoustic transmission into the cabin. 

This is particularly true for flow noise due to the turbulent 
boundary layer over the exterior surfaces. The exciting pressure 
fluctuations are non-acoustic in nature and it is therefore not 
strictly appropriate to use acoustic transmission loss information 
for the skin panels in estimating the transmission, although this 
is often done for simplicity. 
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Transmission loss information for light aircraft and helicopter 
fuselage structures is typically based on extensions of results 
foj. simple panels to the case where the skin is sectioned by 
frames and stringers [10, 11]. Different regions arise based on 
frequency and spacings between the frames or stringers. Radiation 
into the cabin is generally evaluated from expressions for the 
radiation efficiencies of flat plates. Some studies have modeled 
interactions and coupling between resonant modes of the skin 

acoustic resonances of the cabin [9, 12]. These 

studies have sought to identify contributions through individual 
pd.H 0 X s • 

levels for mechanical sources are conveniently measured 
connections between the source and the helicopter struc- 
^ . transfer function relates the measured vibration levels 
to vibration levels on the interior cabin surfaces which are then 
used to estimate the acoustic levels based on a radiation effi- 
ciency for the surface and an acoustic model for the cabin. 

In the case of relatively simple structural elements such as 
gearbox struts, it is possible to analytically quantify their 
vibration transmission characteristics. A difficulty arises in 
analytically describing the gearbox and airframe attachments where 
1 connected. These involve considerably more com- 

plicated geometries. The typical approach for quantifying a 
vibration transfer function is to perform a measurement on the 
actual structure with mechanical shakers, and force and vibration 

S011SOX S • 

A similar experimental approach is also adopted for acoustic 
transmission paths. A speaker is used in place of the source on 
me actual helicopter in order to measure the transfer function in 
terms of a ratio of pressure levels. Where intermediate acoustic 
spaces are present such as the baggage compartment the procedure 
involves more than one acoustic transfer function. 

For gearbox casing radiation acoustic energy reaches the cabin 

through overhead panels and potentially 
transmission into the baggage compartment and then into 
the cabin through the rear passenger compartment bulkhead. For 
the former path a single acoustic transfer function (i.e., noise 
reduction ratio) is required. For the latter, two transfe? 

evaluated first from the cavity around the gearbox 
into the baggage compartment and subsequently with the speaker in 
cabin^^^^^^* compartment from the baggage compartment into the 
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In general, straightforward experimental techniques have been 
applied in generating descriptions of the helicopter cabin noise 
[1/ 4, 5] . Their usage has been successful in provid- 
ing a basis for the design of interior noise control treatments 
resulting in cabin noise levels near ao dBA [8]. As lower levels 
are sought, consistent with the noise environment of commercial 
aircraft, greater detail and sophistication in modeling the noise 
environment and the effects of different treatments will be 
required. 

The remainder of this section discusses the different approaches 
/techniques for source and path identification in greater 
detail . The first^ part deals with the practical engineering or 
experimental techniques . The second part describes more compli- 
cated procedures whose implementation has been greatly facilitated 
in recent years by the development of specialized instrumentation 
and computer systems capable of performing digital signal analysis 
quickly and relatively inexpensively. 


Experimental Techniques 

The goals of diagnostic tests in support of helicopter cabin 
quieting are; 

1. identify and quantify sources of cabin noise. 

2. identify and quantify the paths between those sources 
and the cabin. 

3. describe the resulting in-cabin noise field. 

In the development of treatments to reduce helicopter cabin noise 
it is necessary to create a description of the cabin noise en- 
vironment showing the relation between all of the significant 
excitation sources and the resulting cabin noise level at the 
passenger locations. 

In such a description an excitation source level is an acoustic, 
vibration, aerodynamic pressure, or other excitation acting on the 
helicopter in such a fashion as to cause noise in the cabin, and 
which is not expected to be altered by modifications produced in 
the noise reduction program. The relationship between a source 
excitation level and the noise level at a receiver point is a 
sound transit iss ion path and it is described by a transfer func- 
tion. 
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Its acoustic output. However if th» i a.- ^ atiect 
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a^yroach . - Important characteristics of source oath 


The amplitude spectral transfer function, T ff) relating, tv,« 
excitation at source point n E ff) to the WirilLo relating the 
point k, R (f) is ^ response at receiver 

k , n 


T ff ) = 

E (f) 


n 


( 2 ) 


When N 
series 


sources have N paths 
of transfer functions 


of excitation to a single receiver, a 
may be used to describe the response 


Rk(f) 


f 

n=l 


En(f) 




( 3 ) 


subject to the important criteria: l) that the cower flow aioor 

any of the transmission paths must be only positive in the direo^ 
tion from the source to the receiver 5 \ ZZ airec- 

cross-excitation from one“ soSrle to ' another ^1" r?Lt 
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ducts IS a useful modeling technique. source/path pro- 
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effects on the rotors and fuselage skins. in such an exneriment 

machinery casing vibration would be done to deter- 
ine that this source excitation had changed only a little, if at 

wlJh cabin noise level changes significantly 

speed. It can then be deduced that the dominant source is not 

art proven that machinery and flow sources 

to possible sources, then relating the broad band noise 

tt sources follows. If such exclusion of alternate 

sources does not exist, it is necessary to either do further 

determine the exclusion of those other possible 

excitatio^ a positive relation between the flow 

exciracion and the cabin acoustic levels. 

k“/* methods of identifying acoustic sources 
and paths include flight or ground operation with a single element 
operating (electrically drive an oil cooling fan or a hydraulic 
ii™ stonilrf helicopter is on the ground and all other machinery 
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attJnn^ti'rS,* ’^rv radiating cabin panels with heavy sound 
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Using the simple techniques described here, it is generally 
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Contribution of source/pat h combinations to the cabin n oise 

n of the contributions to cabin noise 
of the various source path combinations is done by a stepwise 

oath^^anH determining the transfer functions for the particular 
paths and combining them with measured, or otherwise known excita- 
tion spectra. Source levels are generally determined from in- 

paJIi,eLrT"'''?rTaV/ appropriate excitation 

functions are generally measured during 
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tative fashion of the in-flight excitation. ^For acoustic e^ci?2- 
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critical. For Vibrational excitation tte dUeluon 
Of force or moment excitation may be critical, as it is also 

Wropriate locations be chosen for the 
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Acoustic excitation . - For the case of acoustic excitation simu- 
lation^ the example of turbine inlet noise provides a good demon- 
stration of some of the considerations and limitations of this 
technique. The turbines are typically located aft of the cabin 
and gearbox. Combustion air is ducted from forward air scoops. 
Paths for inlet noise entering the Ccd>in are: (1) acoustic 

transmission through the duct walls and then through the cabin 
overhead; (2) out of the throat and through the forward overhead 
panels; (3) out of the throat and through the forward cabin 
windows, and (4) through the various other cabin exterior sur- 
faces . 

The turbine inlet noise is characterized by high frequencies (2 to 
10 kHz) which radiate in a highly directive fashion concentrated 
on the axis of the inlet duct. For the acoustic test, excitation 
is generated using several very small, high frequency loudspeakers 
placed well within the intake tubes, giving reasonable expectation 
that the available duct modes would be generally excited. Sound 
measurements are made for the purpose of establishing separate 
transfer functions for inlet noise in the cabin due to trans- 
mission through the overhead panels and through the windows . The 
transfer functions for the two paths are measured independently by 
alternately blocking transmission through the path not being 
measured, using a heavy covering of lead vinyl and fiberglass. 
The transfer function is given as the decibel noise reduction 
between turbine inlet opening and the passenger seats. 

Having the transfer functions, these can be multiplied by the 
measured noise at the turbine inlet source SPL, to determine the 
resulting noise at passenger positions. Separate measurements are 
required to identify the acoustic levels of the turbine inlet 
during operation. The levels are increased by 3 dB accounting for 
two turbine configurations to yield the turbine intake noise at 
the passenger positions. 

Looking back at this crude modeling technique it is appropriate to 
review the following critical approximations and simplifications 
which were made: 

1. Frequency range of excitation limited to 2 to 12 kHz. 

2. Random excitation from miniature loudspeakers presumed 
to represent the spinning turbine compressor excitation. 

3. Effects of forward speed and boundary layer on inlet 
source directivity and transmission loss of panels have 
been ignored. 



4. Incoherent addition, adding 3 dB, used for two turbine* 

toe conteiJuSiHs 

from transmission of sound arriving through different 
panel or window sections. aiiierent 

The frequency range selection is based on fliaht data fnr 
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directional point excitation is likely to result in the rever- 
berant excitation and generalized modal response at points removed 
from the area dirictly adjacent to the excitation. For the 
generalized modal response the distribution of energy for the 
different directions at a given point will remain the same whether 
the excitation is the real or simulated one. In this situation it 
is then possible to characterize the source based on a single 
direction of vibration to be measured during flight and simu- 
lation. 

Local contact deformation at the point of contact between the 
vibration exciter and the structure are likely to occur. For this 
reason it is usually not advisable to measure the source level 
with an accelerometer within the shaker attachment, as is found 
for a standard "impedance head". The location of the accelero- 
meter should be far enough away from the point of contact to avoid 
membrane deformation effects. 

Many of the he^copter components undergo considerable stress 
loads during flight. One place for concern about flight load 
effects is with strut-mounted gearboxes having self-aligning ball 
and socket type connections at each end of the strut. First, if 
there is any play in the bearing it is necessary that the testing 
be done with enough force applied across the joint to maintain 
bearing contact at the appropriate surface. Furthermore, if these 
mounts become rigid against moments under flight conditions it may 
be proper to stop the hinge action by stressing the joint or by 
use of chemical bonding of the joint during testing. Further 
study of these junctions is required before a judgement can be 
made on the necessity of bonding the joint during testing. 


Signal Processing Techniques 

The development of small, high speed computer systems and specia- 
lized instrumentation has resulted in the development and imple- 
mentation of signal processing techniques for studying the dynamic 
behavior of mechanical systems. The procedures involve evalua- 
tions of higher order statistics of the measured response of the 
system in both time and frequency domains . Standard measurements 
simply involve a determination of the system's response level 
relative to an excitation level. 

The following two sections provide an introductory description of 
the techniques and their uses in relation to the helicopter Ceibin 
noise problem. Discussion of techniques with potential for 
application to the helicopter noise problem appear in Appendix A. 




The reader additionally is referred to the following texts for 
descriptions of the theoretical basis for and implementation of 
the signal processing techniques: [13, 14]. 

Source ident ification . - The coherence function provides a quanti- 
tative measure of the relationship between two or more signals 
If one signal represents the output of the system and the others 
are input signals, then the coherence function is an indication of 
the causality between the output and the individual inputs. This 
description requires that the inputs be statistically independent. 

The coherence function, T2(f), is defined by the following for a 
system with an output, y, and n inputs x^: 

|G |2 

-V.? f ns I I indicates 

~G g“ magnitude 

^ ^i’*i complex 

function 



► 

r 

E, 



where the cross-spectral density of the output y and an 
input i' and G are power spectral densities of the 
input and output signals, •^fespectively [14]. 


The coherence function varies between 0 and 1. It equals zero 
when the output is unrelated to the input and 1 when the output is 
due only to the particular input with insignificant contributions 
from other sources . In between values give the fractional amoxmts 
of the output associated with each input. If all inputs are 
identified the coherence function values sum to one in the absence 
of extraneous noise in the measured signals. 


. significant problems with use of the coherence function in 
helicopter cabin diagnostics are related first to too great 
coherence, and second to too little coherence. Gear noise, which 
IS the main source of cabin noise, is generated at the tooth-to- 
tooth points of contact between various pairs of gear wheels. The 
resulting cabin noise arrives through various paths including 
vibration carried across numerous mechanical connections between 
the gearbox and the cabin. 


With inputs defined as the vibration levels at each of the attach- 
ment points of the gearbox and airframe they are likely to be 
highly mutually coherent. The coherence function is in this case 
incapable of distinguishing between the contributions to cabin 
noise due to transmission through individual attachment points. 
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It is capable of distinguishing between the attachment points 
treated collectively as a single source and other statistically 
independent sources. 

The second major problem in the use of coherence techniques occurs 
when there is no dominance of the noise from any single source. 
Flow excitation provides a good example. Because of the spatially 
random nature of the flow disturbances causing the generation of 
noise, the e is little mutual coherence of the vibration response 
of the aircraft window and shin at various individual locations on 
these radiating surfaces. Expected coherence between the cabin 
noise level and measured vibration response at a single location 
will be low, unless there are limited dominating local areas of 
vibration response with major contributions to the cabin noise. 
The situation is one where the distributed response over the 
surface involves relatively small patches of coherent response as 
a result of the spatially random nature of the flow excitation. 

Transmission path identification . - The transmission of energy in 
a complicated system is characterized in general by the effects of 
propagation, attenuation, and reverberation. Propagation refers 
to the transmission of energy without reflection in an acoustic or 
mechanical system. Propagating energy is attenuated in amplitude 
and dissipated as heat. In air the dissipation is due to viscous 
and thermal losses which are generally significant only over long 
distances. In a structure the attenuation is the result of 
internal mechanical damping in the material and accounts also for 
the effect of added damping treatments. 

Complicated acoustic and mechanical systems generally involve a 
great many internal reflections of the propagating energy leading 
to the creation of a reverberant field. The system response is 
characterized by its resonant modes. These are damped by propaga- 
tion losses as well as dissipation at the boundaries of the 
system. For large and complicated systems it is convenient to 
evaluate the overall transmission in terms of individual struc- 
tural subsystems and the connections between them. Statistical 
Energy Analysis is a procedure that has been developed in recent 
years to study the transmission in complicated structural and 
acoustical systems. Its application to vibration transmission in 
a helicopter airframe structure is an important basis for the 
overall modeling of the cabin noise environment in this study. 

At the outset of this discussion it is important to distinguish 
between dispersive and non-dispersive wave propagation. A non- 
dispersive wave is one whose phase speed is ''onstant with fre- 
quency while the phase speed for a dispersive wave varies with 
frequency. Sound in air is non-dispersive as are longitudinal and 
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shear waves in solids. Bending waves in beams and plates are 

^ considered to be non-dispersive if 

the frequency bandwidth of interest is sufficiently naSoi that 
the phase velocity is nearly constant across the band. 

non-dispersive waves is that the 
temporal waveform, moving at the phase velocity does 
not change with time or space. In other ^ords there is nn 

wave^sian?f^^°^^h°^’ origin of the name 'non-dispersive 

or disi,^ri^ initial temporal pulse does not sprLd out 

increasing time. This is because all Sie spec- 

cepLflTv°*'J^d1 propagate at the same speed. Con- 

ceptually, a disturbance in a non-dispersive system can be fol- 

lowed around in time to identify which paths i/ takes to oo 

* . , disturbance to the response point of interest 

^xc‘epf at^tu^5«il:. distu^hance "change 
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INTERIOR NOISE PREDICTION MODEL 
Sources of Cabin Noise 




Acoustic/aerodynamic sources . - Acoustic sources of cabin noise 
include the acoustic radiation from the surface of the gearbox, 
acoustic radiation from the air inlets to the turbines, acoustic 
radiation into the turbine enclosure, excitation of the windows 
and fuselage skin by turbulent boundary layer pressure fluctua- 
tions, sound transmission and radiation from leaks around door 
seals, radiation by the tail and main rotors, and noise generation 
and transmission in the heating and airconditioning system. The 
acoustic radiation from the gearbox is the most significant source 
listed above. However, in a fully- treated cabin the other 
sources, in particular TBL excitation of the windows and doors 
leaks, can be major noise sources. 

One of these sources, airborne gearbox noise, is caused by mechan- 
ical 9®sr mesh forces which generate vibration that is transmitted 
through shafts and bearings to the gearbox case. The vibrations 

9^ gearbox case radiate noise into the acoustic space enclos- 
ing the gearbox. 

The acoustic power radiated by the gearbox at a particular fre- 
^ency is related to its vibration by the radiation efficiency, 
the characteristic impedance of the acoustic space enclosing the 
gearbox, the surface area of the gearbox, and the mean-square 
velocity normal to the surface of the gearbox averaged over the 
gearbox surface (see Equation (1)). 

The radiation efficiency for flat panels is a function of fre- 
quency relative to the critical frequency, f at which the 
bending wavespeed in the panel equals the speed of sound in air. 
The critical frequency is dependent on the material properties of 
the panel and its thickness, h. For aluminum, steel or glass 
panels in air, the critical frequency is given by: 


= 12 700/h 


Below the critical frequency 


a 


rad 



( 6 ) 



where p = perimeter of the structure 

= acoustic wavelength at critical frequency = c /f 
P = panel edge boundary condition factor 
P “ 1 simple edge support 
P ^ 2 clamped edge support 

^al frequency the radiation efficiency is generally 

critical frequency due the radiation below th^ 
For each rib^f lenath L the perimeter. 

of ''the'''tib’'‘’® D^Micall^^'"'^ ■7°”' 

blade passage freguencies uhiio k .aw rotation rates and 

primarily f?L cStioL broadband random noise results 

™easu«d"i[b?a?fon"1“elTi;e:- 

r?,f eo:^«‘^Ve?'?lr" the“turbii.r ation'lSli^:": 

measurement of sound power in \ free" f1 el 2°^®® • ® 

lently a measurement ”01 reve"rb^/a"nV sou'nd I^^^^S^Tev^el 
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aiffxcult co separate turbine inlet noise from airborne 
gearbox noise or airborne turbine casing noise. airnorne 

cSbin^a?! transmitted to the 

A ci rrri^ direct, and this suggests that the turbine noise is 
hv^ a significant acoustic source. This suggestion is supported 
y e noise control program for the Agusta A-109 in this 

Kg^mcaf/""but"\^^^^^^ wai" found "J^be 

igniiicant, but the airborne turbine inlet noise and i-ht-k-; 
casing noise were found to be insignificant turbine 

A third source of airborne noise to be considered is the turbulent- 
boundary layer. At the higher flight speeds toe pressure fluctua 
tions in the turbulent boundary layer ar^ a soGtoe tf bMaSt.nd 

;?Srs l?/"oute°riAin ^■'/^toe4^%i^f%.^^“i“=. 

indirlct^JSdLtioS^f?o^^t^^°” Sindo^s and 

aJI^ifed!"”^ compared to the window radiation unless double windows 

speeds the pressure fluctuations in a fully developed 
ampiricahv r?f ^ ^e priL^tJd 


rms 


= 0.006 q 


(7) 


where is the free stream dynamic pressure, 
% = pUj/2 


( 8 ) 


Where p is the air density and is free-stream velocity. 




( 9 ) 
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where 6* is the displacement thickness of the boundary layer 
Below • • - . - . J ' 


“pk spectrum is not a strong function of frequency. 


However, above lu . it falls off approximately 6 dB per octave 
(u) 2 dependence),*^ 

The boundary layer displacement thickness depends on the geometry 
of the flowfield. Given an exact geometry it is possible to 
calculate the displacement thickness over the surface of the 
airframe. Because of the relationship of surface drag to this 
parameter the calculation is usually carried out as part of the 
design. However, as an approximation for the purposes of 
cabin noise prediction, the displacement thickness can be de- 
termined from the equation 

6* = 0.0016 D (10) 

where D is the distance from the leading edge. 

Tile effectiveness with which the TBL fluctuating pressure field 
e.icites the window or fuselage panels depends on the correlation 
lengths of the field. The spectrum of the modal force, S^(u>) is 
related to the spectrum of the TBL pressure, Sp(iu) by the elation 

Sj(u>) = A2 Sp(w) (2j) 

where j is the joint acceptance and A is the window or panel area. 
The joint acceptance is determined by the relative size of the 
correlation lengths of the TBL pressure fluctuations and the panel 
area. It is given by the equation [16]. 


j2 = 0 ^ 1 

71 A '•Ui ^ 


The panel response is found by computing the response of a typical 
panel mode and multiplying that response by the number of modes 
per unit frequency. 

The response spectral density averaged over a band of frequency 
that includes several modes is given by 


= f J 


i“ Sp(u)) 


0 





radiation i^to \he c^ir cVn ca?du?alel“fro» J”' 
spectrum and the radiation efficiency. ^ vibration 
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^^io n somrcGfi • 

^e most important contributors*°to modern helicopters are 

®*ivironment. The dominant vihra^ • cabin interior 

which IS the primary focus of gearbox 

factors contribute to this asses^^nf ® ^ number of 

of mechanical power avail^le? the large amount 

mechanical power feeds throuah from that the 

considerably lower speed rotor turbines to the 

nature of the helicopter desi^ 'conflgj^atio'n 

^ imperfect meshing of the aear -i 

f H° ® excitation of the gear J^® steady loads 

at fundamental and harmonic tram and gearbox casing 

mash and shaft rates. TSe vib^aS^^n ^®^ “ '*® Particular geal 
through the bearings supporting the ‘"""smitted to the casing 

study IS not concerned w'to tte detailed P'^«“nt 

the generation of dynamic gear m^ach mechanisms involved in 

mission of vibration to tie ge« ca^^^^^ trans- 

defined in*terms^ii%hJ\ibrftion^^ gearbox is 

with the airframe. For strut-roountad*^^ ^^® ®^^®chment points 
considered to be a ^arf oJ ?an bl 

g sources is one which 

to which the source is attached structure 

mann®r allows for modifications^of the attach”? ^?® this 

the necessity of changing the source without 

of importance for the Dresen-h ^^P^^sentation, An example 

isolation mounts. Independentlv ^^® of gearbox 

for evaluating the effect of th» ®P^®®®tting the source allows 
to changes in the sourL levels for ™e geartox’; due 

tional sources in'^^three'liirections ”*’'** transla- 

presented either by blocked ^ 1 source level is re- 

The blocked force ^eve? !t tte forc^® ^'^®® levels, 

point to identically constrain i-h<» ^®g^ired at the attachment 
vibration level i«s I- t motion to zero. Tha 

disconnected and the forc^ is zero'!^^ result if the attachment is 
to^ t^e gear .esh forces/„otioL^Vt"V*''gea^" ?e%"t'h“i^ 






wriGiTG ^ V f* ' 

forces, ‘respictively,' in\h«e“LectiMnt^tir®i/*i°°‘''^*= " 

points wfthin the gear/Sx wh«e ?he®i'i‘'h''‘r® « th 

The complex mobilities (ratio of velocitv*' generated 
^12 ... Characterize the dvnamir applied force) Yt , 

(gears, shafts, bearings, LrcSsin^) ° gearbox structur 

« "Obilit. 

the attachment forces f ^ , f , , f, tre ze?n velocities when 

general not equal to zero Thi g®ar mes] 

free velocities at the attacl^nt point defines a set oj 

I T7 free 1 


free 

free 


I Y,4 • . . Yi^ 

j Y24 • • • Y2 j 

( Y34 • • • Y, 


( 15 ) 


is possible , 

point in terms of the free velocities- attachmer 




which can be written as: 


[Y] ? + ^ 


free 


Tterattached structure can be characterized by an impedance 


I = f 


= - <z 


‘11 Zi2 Zi3 

■2 1 Z2 2 Z2 3 


^^2 } = -[Z] V 


( 17 ) 


:f"\he^a\4Thed"rtlucii«^‘/n Vrm's '^*'’avior 

attachment point. The minis forces/velocities at the 

where velocities are positive into convention 

Velocities Vi, Vo Vo are definoH structure, 

therefore the veloci\ies in^i ® positive into the gearbox and 
these. Equations negatives of 


V^^®® = (1 + fZ][Y]) V 


(18) 


atiachld®llructu?r‘'®rK impedance characteris- 

1) the attachment point structural^^impe^anLs^iTd^^J ^ 
behavior locking into the source dynamic 

described by the mobility matrix fYl point as 

structure in the form of an icr^^ =.+-• ' Modifications to the 
changes in the impedance mftrlx result in 

-d ^ attachment point 

may be generated for oth2r S.^"'^lar representations 

angular velocities and for oth« s^i'rces *S>°f vVbrato^ moments and 

possible to calculate free vein/J??,, i' ,^^at case it may be 
levels based on Eouation flfii an^^i-K ™®asured attached 

the source and dynamic characteristics of 

velocAy levefs can °"'® calculated, the free 

structural modifications *** evaluate the effects of 
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The solution to the overall 

of information relating the velociti”e/^^^ 4 -K®® additional piece 
the cabin interior noile levets L attachment points to 

a set of velocity transfer functions from^^h^ to define 

on panels that radiate intn ®^^*c^ent points to 
vibration levels are combined with a The panel 

mg the effects of noise tre^iSini-= ^ efficiency, includ- 

acoustic field to comolete ® model of the cabin 

velocity transfer fSnSo^s for the 

field are described in following section^^.^^® acoustic 

gearbox vibrVuo^n transmission^ concern in modeling 

and transmission paths. The internal fo coherent sources 

the gearbox are a source of coherent mesh mechanisms within 

attachment points Also the vibration at the different 

ment point to f vibraUng p^nt ^ 

passenger position in the cabin a«ustic^ffela®o®^ surface or to a 
but coherent transmission oaths t?i^ ^ occurs by different 
ence effects in algebraicallv summing ignore coher- 

contributions as though they were statisticafly 

''®The7rerve!ocit™s Egua- 

event that the appropriat^^ mobiT tii r ^^^°^®rent in the unlikely 
so that the inte^liai fo^Lr^^^ are zero 

velocity each. This is extr#anf*ai t/ i only one free 

ternal structure of the qearbox considering the in- 

expected to display a high degree o? cohl?en?e . 

diagonal^ te^rms the mobimy^and^^^ non-zero, off- 

tions (14) and (17), respectiveTv matrices from Equa- 

tached velocities v, Vo^ v wli ^^^^3 in coherent at- 

points on radiating cabin' panel surfacet^ functions to 

ferent coherent paths. Th^ ovel^ll eff^^i ^ involve dif- 

coherence is mitigated by the fact th^^ i ^ source and path 

levels within the cabin are of inSrLt ^ average response 

djlcrtb^^d^ has been 

than the gearbox, such as the hvdrainff^ sources other 

mounted directly to the airfrrimJ^ lines. These lines are 

between the pump, which is monnTo/^ approximately 0.5m intervals 
for the -in"anl'taA^?oti.^rsTnrt"he°US!nne"^?*- 



Applying the above procedure for th^» hv<^r-ai,i -i «« 4. 

attactaen\"s^Te fdeSSf “d *:otioTs' "a? 

a???sf a'^d FiPPJ^^ 

=„t^ora\\a?S 


sysLST^^LIeMr^J” ItUTtS. ^oh «%4“™'ed°U^the?‘‘^lSJ 

represent an overall powe?^ balaLe^ 

vibratory or acousi-ir'^mrt/^^ao each subsystem. Resonant 

energy JlthYn «?h“;Ssya ttm “* >^«Ponsrble for the storage of 

s“%rSVZ;‘£s ’£4 ““ •• w«„“= 

?S2^r*iTtnVe\“ 

f°^de?all ty"°^1,y, batfr“?h^*oV*EISL^‘l?A-^1= 

?og"etier *‘^°"' modes 

t^^;odr-S“^^ 

finrte element or modal ana?ysis S the°"aianfbf^^^^^^ 

anSTo^sf al^e 1cVi“^ed“ 

structures such as a helicopter. model for complex 

=S;iHHES“ “ 

:::r„i- .r ““““•“ 

allowed to have different energies different groups are 

fJlm^tUrfield of “roo^a^oSst i“/s "°1T\ 

acoustic space can have thnnaena'. concert hall or other 

range of interest. To stSd^r/h f crprexTs^emn*cou^J??c“S 
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use the concept of a diffuse field. In such a sound field acous- 

point comes from all directions, with ‘»nerav 
from different directions being uncorrelated. This wave descri?^ 
tion IS equivalent to saying that there are a largr nurler of 

is^^use/^?<!? the same energy. Again, equipartition of energy 

IS used to group modes together and to study complex dynamical 
tion^"'^ without having to identify each individual mode ojTvibra- 

modeling procedure, the success in using SEA depends 
on the skill with which the model is set up. Often it is advanta- 

0 ?°^ shnolTff^H^ predictions with data from a laboratory study 

approximates in some general way 
the actual structure being considered. The goal of the laboratory 

1 ° obtain data that can be used directly to predict 
suLort ^®sponse being considered. Nor is the ^al to 

^ of SEA. The true purpose of the 

modeP th7t has\een lit validity of the particulL 

anui ^ % 5 the model can be 

applied to the complex structure being studied. 

v” using SEA depends in large part on grouping the 

samI eieJav^^^The^fn?! within a subsystem has the 

same energy. The following principles are used: 

J • L subsystem should have resonance frequen- 

cies in the same frequency band; 

Within a subsystem should be identified with a 
particular structural subsection of the system, i.e a 
plate, beam, or shell segment or an acoustic space; 

® particular structural subsection should 
consist of the same type of vibratory motion, i.e., bending 
compressional deformation, shear deformation, or torsion'. 

In highly complex structural subsections the modes may involve 

vibratory motion. in this case a!? modes oFthe 
subsection having resonances within a frequency band can be 

i“ Le group For 
^ section of a beam has three groups ofmodes- 

mo?for longitudinal motion, those involvi^^ torsional 

motion, and those involving transverse motions. The transverse 
^^urther divided into two groups corresponding to the 

Its length, SEA allows two approaches. in one approach each 
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straight section of the beam is considered to be a separate 
subsection, and the four groups of modes in each beam section are 
allowed to be coupled with each other and with modes of other beam 
sections. A second, far simpler approach, is to treat the entire 
beam as one subsection and group all modes into a single group. 

u ^ possible because the bends will result in normal modes for 
the beam that involve all types of motion - longitudinal, tor- 
sional, and transverse - 

Because these modes consist of all types of motion they tend to be 
similarly excited by external sources, similarly damped, and 
similarly coupled to modes in other groups. Under these condi- 
tions the modes will tend to have egual energy. 

An SEA model consists of a block diagram where the individual 
blocks represent groups of similar modes. A sample block diagram 
and the system it models are shown in Figure 18. The system con- 
sists of an acoustic space, a panel which is one of the bounding 
surfaces of the space, and a beam member which is attached to the 
• The configuration is representative of subsections from a 
helicopter airframe, but is by no means indicative of the full 
size model required for the complete airframe. 

The groups of modes or blocks in the model include the acoustic 
resonances of the space, the bending modes of the panel, and 
bending modes of the beam. A more detailed descrip- 
tion of the modeling of subsystems of this type is presented in 
subsequent sections of the report. Lines connecting the blocks 
represent power exchange between the different coupled groups of 
modes. Additional lines represent the power dissipated within 
each group of modes and the power input from external sources. 

The input power in this case is due to other structures which are 
attached to the beam. The dissipated power is due to internal 
mechanical damping for the structural subsections and wall absorp- 
tion for the acoustic space. The transmitted power is the result 
of mechanical connections between structural subsections and 
acoustic radiation into the acoustic space from the vibratinq 
structural subsections . 

The identification of power transmission paths is initially 
straightforward. Mode groups from adjacent structural subsections 
are almost always coupled, although cases can occur where this is 
c^se . For example, in the case of two beams joined at 
right angles, the longitudinal modes of one beam do not couple to 
the longitudinal modes of the second beam. However, this is the 
exception rather than the rule since the longitudinal modes will 
be coupled for any other angle of attachment. 
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modes ^of*^struo\?ral°*ubserti^^^^ S‘®arr^°'’l between 

another through a third interJ?nin!i o k connected to one 

of frequencies it is"" por^Srr for ^ ^and 

structure to be coupled to the resonant modes of one 

ture through non-resonant modes of T modes of a second struc- 
The most common case in which this n^t^, intervening structure, 
important is that of two acous^?^%L^ transmission is 

shell structure. The modes of separated by a plate or 

coupled to the resonant moSs of th. directly 

to the resonant modes of the other acSnQi-^^^ indirectly coupled 
resonant mass-law controlled mode? of through non- 
path IS often more important than indirect 

3 ^ndi 1 1 ons °Le ^derived” by^lq?a??nV the t^ tor steady-state 

a group of modes with the sSS of th^^J^m to 

and the time-average cower transmit power dissipated 
For example, powef balance f"r thl of modes, 

acoustic space modes in Figure ?e%an be”Srittira^s 


diss 


_ rj(3,4)_ 

trans trans \rans~ ° 


(19) 


where 10 +-Vi 

group bhe^net Power transmitted from mode 

and W<f ie the time-iverage Vw-e^ ^di'^st''a\ed" 
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balance equations depends on the basic SEA 
A . ^^®®"®verage power exchange between two groups of 
similar modes is proportional to the difference between^ the 

transmitt^i?^^ group. The time-average power 

transmitted from group i to group j can be written: 


N N f _ _j .tot , 

trans l u ». J 


coupling coefficient between a mode 
^ ^ ^ number of modes in 

in qIouv i^^^Tt^Wp time-average energy of all modes 

ticf] rnnipt^+-^!?K average coupling coefficients is a statis- 

tical concept that is an integral part of SEA. 

^ includes the 

^ ^ resonance frequencies be within a given 

E^ftifn^'oor ^iven this requirement, the coupling tei™ in 
Equation (20) is almost always replaced by a term containing a 
coupling loss factor and a band center frequency so that: 


= ujn - . N. 
trans 'ij i 


E. 

r^>tot 


factor center frequency and n.. is the coupling loss 

the hanH ^ time-average power, iV now the power within 

the modi^lSInts , N. and N. , include 
only modes with resonance frequencies within ^e band^ Aw and the 

the^band^^^A^I' tQt, include only those modes within 

ar^'iffowed to be functions of frequency 
frequen^cies^ dynamics of the system can be studied at different 

'usually taken in evaluating the transmitted 

modal Hpn^^ vf ^1' evaluated in terms of a 

odal density, which is the average number of modes per unit 

frequency in an ensemble of systems in which the resonance fre- 
quencies are randomly distributed. The concept of a^ ensemble of 

sistint with distributed resonance frequencies is con- 

<R>^ average mode to mode coupling factor, 

fimAi-i/^r^ modal density is a slowly varying 

be Written count in the band Aw can 
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= n^(iu) 4uj 


Where n.(w) is the modal densitv at th^ k ^ 

in casei where this is not tr/e, EguatiorljIT^i/ r!p?ac27b;.. 


~ / du) n. (uj) 

Au) ^ 


install o?"|',u«1o:‘‘%rT3“n"ot°'’nf 

balance can be written: ’ * necessary. Thus, the power 


0 


= lun 
trans • 


• . n . 
ij 1 


E. 

r i.tot 


j "" m\nu^ "’th^ transmitted from 

with Equation (24) to obtain the^relultf ^ ^ combined 


Oil n = Q. . n. 

^ (25) 

loss factor, since i\ is^sometim^ calculating the coupling 

coupling loss factors, n . Z I calculate one of thf 

be related to tL^totSl^eneJgy^in the^^ro^ modes can also 

definition of a damping loss fac^ or ^ i°ss factor. By 

power in a band of frequencies ?s g7v'en by: dissipated 


y(i) 

diss 


WO. E . ^ 

1 i,tot 


re«!ttA il 


Equation (26) is usually 


y(i) 

diss 


WT| n. 
1 1 


iLitot 


0 


so that it is in the same form as Equation (24). 

balance egua- 

tions, a set of linear algebraic equations can be obtained vnr 
the system in Figure 18 these become: obtained. For 

Subsystem 1: Bending modes of the beam member: 

^ - ”.n,3 ^ -n,n,, ^ = «.<■;» (28) 


Subsystem 2: Torsional modes of the beam member: 

E, 




' 3. tot 
'23 n. 


" 2^24 


4. tot 


n, 


= 0 


(29) 


Subsystem 3 : Bending modes of the panel 


« h ,tot _ ^ ^2, tot , 

3^31 - "3^32 n, ”3^ V'^31^'132^^34^ “7 


~ n n 

3 3^34 ~ ^ 


(30) 


Subsystem 4: Acoustic space modes; 


_ _ l.tot 

-t ''aH 


1 


"4"43 = “ 


latOt _ ^ “3. tot 

4' '42 n. 


(31) 


Figure 16, which likely involves many hun- 
fmir^ -i^d Structural and acoustical modes, requires only 
dSnLii subsystem energies in the SEA model of its 

and are re^dTiv^* above equations are linear and algebraic 
the readily solved by standard computational procedures once 
the SEA parameters have been determined for each subsystem in the 

interest. The SEA parameters include the mSe 
density and coupling and damping loss factors. 
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model'^with many of systems can require a 

power balance Md ?h2 ?elatfonshir, time-average 

variables continue to be valiV^ lolution'^o'V *“''9^ 

equations can be tedious Hoi^etrsmv- power balance 

inversion routines with a small dioit numerical matrix 

be easily obtained e"ven^L lar^ge " ys^'Ts'*.'""' " 

aS^a„\aTe "ol'^si^g S%^“moTl“ ^ “sjor 

and analytical expression haJ»' has been set up 

power inputs, loss factors OymOf," “Stained for the required 

quickly obtained for a variety of di f feren? can be 

ariety of different design configurations. 

iill^^L^dev'Jlopedf^w^^^^^ SEA modeling 

involving typical mechanical structures "Lralou^stYc''^^^^^^^^ 

^stm ^gr^e^^ des~cr?by^b^ryerfei°^ ^ ^ 

a?^0t?aightfo0,ar?lflx"tenl:5 the" particular “system'Yn V^|u“re 1 
for the geneOT s«Lriro0a0O’1f„®'^\ Each eqLtion 

vibratory power iroJt to a / balancing the time-average 

resonance frequencies within a hanH^f vibration having 

average power transmitted to IthVl J'"’ ”"bh toe sum of the time- 
average p'ower dissipated wiOinOhe StT “4.'" It® 
variables are given by the following^xplessionsf 


Input Power Spectrum 
Dissipated Power 


Transmitted Power 
to Mode Group j 


ujq.E. ^ 

1 i,tot 


I 


tor|rot^*from lu ««ces'e««la?^to‘’?^rs i" 

density for the arouo u, ic system, n. is the modal 

dissipation loss Sr Vr%Jl®mo^®e"s ^ 

average energy of all modes in J/tot^^® totar time- 

g?21p"rtl"grlir^"® time-avera^ po«.r^^^n%nitltTtrol 

«ttointoVTto?al “®m mod\ to lie" eguiult ®"“®” 
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■,<n, . ■ , ) ,JJa, 


■' j, 


ij'‘ n. ■' ti) \ 

J i 


ship can be written: ' following reciprocity relation- 

i ~ n.n. . 

^ IJ 1 *1X 

( 33 ) 

ing^thl^eneTgy''^|er*mode in®tach ® sy™»etric matrix relat- 

iipear set o/^egSatioTs caif thif be"Tftte°n a^^ Power.^^'lJe 


•^inii 

"^202 1 

"^sOsi 

— 

— 

'"2021 

02022 

“^3032 

— 

— 

■"3031 

*^3032 

^3033 




■n n 

m 'ini 


El. tot 


•n n 

. m 


n ^ 
m nun 


m, tot 
n 


where the total loss factor is defined as 

“ '’i * jh "ii 


0 



7 ^ 




ESv3‘13> 

s«?iEr --- 

with a minimum of computational effort. ^ energies 

^®="i'^i''9*'thl ftn®d^ntLx 'Approach' 

and“coi"la'?g l«s®4cto?I . 

a.w aspects or SEA. A subsequent section describs »<5 i-h® nd 
sfythese expressions for the particular SEA model of the Sikorsky 

ggi y y a-a/:; ;:z:r;.g;r a s:: 

-s> "“n&4" z “-s: 

which a model Xr^^cture is eSted“hv® experiment in 

::r." 

js*^piikh«'£f IE r ^ 

components and thereby usrthranalnic” iSd?™Snt* 
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The analytical expressions are given in terms of a mode density 
n(w), or average number of modes per unit freguencv Resonan4 
modes of structures occur at discrete frequencies? and! ‘ therefSJe 

will?xhfbit stSTdi^ ® particular frequency band 

in the band t>^ additional modes are included 

analytical expressions for the mode density are 
smoothed estimates of the number of modes in the band. 


SrdiTg^\o?Sffofloilng?'^^ "i<“> 


u) + Aiu/2 

N. = ° f 

1 


u)^- Auu/2 


<iui n^(ui) 


( 36 ) 


In most cases the modal density is a fairly smooth function of 
frequency relative to the bandwidth of interest so that 


N . = n . (u) )Au» 
1 1 o 


( 37 ) 


expression for the mode density is derived for 
Ideal structures having ideal boundary conditions and precisely 
accounts for the location and distribution of the re sonLi- modes 

Real structures are characterized by non- 

structurr??ofL random non-uniformities L the 

tructure itself. These act to shift the resonance frequencies in 
comparison with those of the ideal structure. r.equencies in 


uilikelv ?rf I!?? density expression for the ideal structure is 
unlikely to give the correct number of modes for the actual 

acJu^Sr'fof -timate. That estLa?;' is'„o?e 

accurate for broadband excitation at high frequencies where the 

?r’narrlJ"°h1f„^ predicted to be large.® At lo®wer freVencies^l 
barrow band or pure tone excitation the uncertainties in the 
mode density and SEA estimates of the system response are greater 


st?uSti?ei dfnenHin^ as one, two, or three-dimensional 

?f inte?est?p1«?i^!i T ® wavelength of sound at the frequency 
rCi i t to cross-sectional dimensions of the space^ 

Where the maximum cross-sectional 
dimension is less than the wavelength of sound are modeled as 

Spaces. The mode density depends only on the 

length of the space, £, and the speed of sound, c 





n 1 (i") = i/ (nc ) 


( 38 ) 


Thin, flat volumes where the thickness is *-k 

wavelength are modeled as acoustic 

density is given Jy, two-dimensional spaces. The mode 


= Au)/(2nc 2) + P/(2nc ) 

' U 


(39) 


area. The^ mode Snsity for a^three^di ^ perimeter around the 
given by aensicy tor a three-dimensional acoustic space is 


Ha 3(10) = U)2v/(2n2c 3) + uS/(8nc 2) + p/(i6nc ) 

O 


(40) 


length for the^spac^^' ^ surface area, and P the total edge 

subIJsLms^®LclLi%^°™o^^^ importance for structural 

motions. The propagation wavesoeeds' - and flexural 

longitudinal and torsional different. For 

stant with frequency a^s the wavespeeds are con- 

Such systems are referred t? as hfs the propagation of sound, 
dimensional systems the non-dispersive. For one- 

with the sound speed slt^ ecLai P^^opagation applies 

longitudinal or to?sfonarmotTn! appropriate value for 


= £/(nc) 


where for longitudinal motion in rods, 

c = 

E = Young's modulus 
p = density 


(41) 


(42) 


Whereas for torsional motion of straight bars, 
c = VGJ7Tpip) 


(43) 
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G = shear stiffness 

~ torsional moment of rigidity 

I = polar moment of inertia of the 
cross-section 

ttl polar 

the shear wavespeed for the liat&'ial ’"“^sional wavespeed equals 

The moment of rigidity, J, is discussed in references on the 
is gi?en by '1®)' ® formula for approximating J 


^ 

1 + 1 

I T 


th^x ’sn^r inertia of the cross-section about 

action ?p.^*a!s ® ‘“e cross- 

^ (45) 

For a square rectangular bar of side A the moment of rigidity is: 

J = ^ I = T - ^ 

6 ’ X y ■ 12 (46) 


For a narrow rectangular bar J becomes 


J = ^3 

3 


b - long dimension 
a = short dimension 


0 
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moment . of \igidl?[^f VcoL^Png ^elaUve^^to'^'^he 

no-warping condition becomes: ^ lacive to tne 


wa rping 
J 

no“wa rping 


4a ^ 
b2 


« 1 for I >> I or b>a 

X y 


(48) 


For aircraft or helicopter structures vibration transmission 

iSr /r\mes members is often important. 

moLnts o? ^ 5^®“ angles, channels, or I-beams. The 

thr exoreLio^ for 5hf structures can be evaluated using 

narrow rectangular bar. The frame is 
rectangular sections, and the moment of rigidity is 
approximated as the sum of the moments of rigidity for the indivi- 
dual rectangular sections. Alternatively the moment of riaiditv 
can be evaluated from Eq. (44) directly ihere anT 1° applj fo^ 


the full cross-section. 


For bending or flexural 
wavespeed is a function 


wave propagation in beams 
of frequency: 


and plates 


the 


c, = C„Ui 
b ’ r £ 


= radius of gyration (49) 

of cross-section 

^£ ~ longitudinal wave 
speed in material 


For a narrow beam: 
= .JeTp 


(50) 
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For a plate: 


^s. = Vk/Ip(i-o^) j 


= h2/12 


a = Poisson's ratio 
E = Young's modulus 
p = density 


h = height of cross- 
section 


For a one-dimensional structure of length £: 

n (w) = ~ (i-iii 

s,l Ttc c dw^ (53) 


and for a two-dimensional structure of area A: 


plates these yield, along with Eqs. 
following equations: ^ 

For a beam 


(50) - (54) the 


n. (u)) = — 
beam 2nc. 


and for a plate 


n T (lu) = — - — r 
plate 4 tik ci 

r £ 


0 


a.slipa«a*wimra ' subsystem s?orel’'en\rw* 


= lun. E. 


diss 


i i,tot 


( 57 ) 


eS«;y> 0 ^^^" iir^inc1Sd%°\r dh.fSSt '^^ the stored 

Within the bulk materials of ±h>^ o^,k aisSrp^tion mechanisms 

at the structuraTco^nec^^^^^^^ 

* iT3,di ^ tod by vibr t' f*oir som© problems is the acous- 

acoounterfor ifthe SE^mSdeJ”® "'^™=^“tes which is otherwise not 

le^el!^^ For*’^m^taTs° tht^ at the molecular 

5Srue“’r li-caJS-aS !r “ 

tmr ■ so??lr‘^%LlS:e1?c"^^ll“^ valilfin^^thr^^re*- 

ranging from .1 to .5. have higher loss factors 

trlaSfs' "/uc”\«%tae®nts®oL\%fe??'’°'’t“'^® “ 

of the panel to which thev increase the damping 

filtorl “c^uX X^^thfctne^-'^^^^ie^^i 

miric “Itir^afs to™fr ^amoLr^^a^^™ Pt°tesses within elasto- 
tions of both temperature and elastic properties are func- 

treatments typTcaUv mSrtai^ ^^^^^^^turers of add-on 

frequency and temperature r^noAd particular 

determining composite loss f^toVs for provided for 

to different thfckness panels treatments when added 

damping values ^a^re ^gene?a\^lV^ l°arger^ tharf structures the 

itself as a result of damoina H?i^fo ''^lues for the metal 

fixtures, etc. Dampina values for connections, attached 

to have a broad range of values near .01 "^^^ctures are known 
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Two theories are proposed to exni;,-;^ 

riveted structures . One theo?v k greater damping for 

and involves plastic deformaS of slip Wl 

joined surfaces. a aas ^ ^°«tact asperities of the 

source of the dissipation^ for Proposed as the 

losses are generated as gas in the °aan^ [20). Viscous 

surfaces is pumped back and forth hv^i-h , 'the contacting 

faces. The air pumping mechanism "^^e sur? 

damping estimates consistent with provide 

aerospace structures. measured values for riveted 

Pp\rt^\“L“"irtabYisSr/",3>a^^^ ‘ha structure is 

damping treatments are often a fir^r p. levels. Add-on 

reduced vibration transmission aiid nois^^^'^Tlfi' 

that exists m estimating dampina loss precision 

problem for all vibration\naT^s?l tl“c?5.i,"u%T[Ll^u^d/„,1|^!““-'‘ 

absorption. *°Ai“‘^so'^Uo^f“inv^l^^^^^^ result of wall or air 

energy mto heat as sound pronaalt^.. "=°^^®^sion of acoustic 
acoustic spaces in the audio freauencv%^^''°''^ smaller 

typically significant in compariLn Sit^ wIll^i^sSp^i^^"^^ 

on reflJction^Sf^the^walls^of^the^ in energy of the sound wave 
that include the resorant^odes of For SEA models 

adjacent acoustic space on i-hra r^+-f the wall structure or of an 
the energy transmitted to these mo^es^^"^^ of the wall structure, 
pation in the source space bif i<= “odes is not included as dissi- 
factor. wall absorptfon accou^L ® coupling loss 

at the wall surface^ as a result ofl^cf^ ^^® dissipated 

resistance of the wall. Abs<fbira 

carpets or fiberglass blanhets conVi^\e%^o^“^ L^^r^^Ton" 

are nof in<fuded modef^tL^^^ adjacent acoustic spaces 
wall structure and tS?ougS Tt to Ht ^/^^^“‘i^ted into the 

then accounted for as dissipation in the sour^ sfacf 

absorp?!on"coe«ic"en??\-“^ the 


c S 
_ o 

n - ^TTTT — a 


8;tVf 


( 58 ) 
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where 


- speed of sound 


S - total wall surface area in the space 
V = volume of the space 
f = frequency Hz 


The weighted absorption coefficient 
weighted average of the absoSiSn 
s^J^f^ces in the space: 


is determined from an area 
coefficients for individual 


a 


o ^ S.a. 
S i 1 X 


(59) 


where 


^i ~ area of individual surface 


Of 


absorption coefficient for individual surface 


sSrbTng^ ^P^'^ticularly those involving ab- 

significantly graater^an Se en^ 

to an adjoining acoustic sm« or through the wall 

wall structure. Reverberation i- ■ resonant modes of the 

accurate estimate of th? energy Provide an 

measured damping loss factor walls. The 

be used directly to represert th#» can, therefore, 

equation for the acous'?ic%Tce."'® <»i=^iP*tion' in the seA 

e“lSating“'"S''?f£e°ct o? an* ®"® interested in 

acoustic space. Data on ^existina^trea^^^®?^”^ treatment in an 
form of an absorption coeff^cieV generally in the 

damping loss factor for the space wi^h^-h ^® 

ment is evaluated from Egs (58) and absorbing wall treat- 

is characterized by the ahcsonni-l^^ treated area, s. , 

absorbing material, a . The untreated add^d 

by the average absorption for the 


( 60 ) 


+S^Cf. 


“t = 


s + 

ut t 


eS® "*‘® treatment is then evaluated from 

io^siA"*^ ■ ’’’’® coupling loss factor (CLF) is unique 

K ^ * characterizes the dynamic behavior at the iunctiors 

between connected structural and acoustical s >.systLt Tf 

* Specifically the coupling loss factor determj s^c 3 
^lp%h across the junction relative to the storoo energy 

of the source subsystem. The calculation of coupling loss factoS 
can be the most complicated analytical step in SEA, and there- 

experience and understanding of the' dynamic 
behavior of mechanical and acoustical systems. aynamic 

important in determining the distri- 
svsteL ^ i energy levels amongst the different sub- 

4-^ system with uniformly distributed damping levels 

thi energy levels will be greater near the source for 

the case where the subsystems are weakly coupled, that s when the 
coupling loss factors are small. For'^large CLF -s Moratory 
energy is distributed more uniformly throughout the structure^ 

to^rGDf^ee''Trf'' analytically evaluating the CLF is 

ihill^maintainnne’^1-^^^ Structures with semi infinite structures 
^ junction geometry. The power trans- 
mitted across the junction is evaluated for travelling wave 

poSer^t? the^w subsystem. The ratio of the transmitted 

??a:Lissiol! 00“^^^^^^ ^“veiling waves defines a 


_ trans 

n. 

me 


( 61 ) 


o? systems such as frames which deform in bending 
or torsion the transmission coefficient is related to inout 
impedances for the connected subsystems according to: 
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( 62 ) 


4Re(Zj) Re(Z 2 ) 
^,2 = II ZJ2 


where 


^2 = input impedance of source, 
receiving subsystem 

I 2^ = sum of input impedances at the 

junction for all attached 
subsystems including source and 
receiving subsystems 


The impedances are either for force or moment inputs depending on 
the geometry and the type of motion allowed at the junction. The 
above expression applies when the transmission is due only to a 
single type of junction motion, either angular rotation or linear 
translation. 

Accounting for the difference between the transmission coefficient 
(where the transmitted power is proportional to an incident power) 
and the coupling loss factor (where it is proportional to the 
stored energy) results in the following expression: 

c 

^ 1,2 ~ 2 w £ ^ 1,2 (^ 3 ) 


Cg - group velocity for the particular type of motion 
Cg = 2Cj^ for bending waves 

Cg = or for torsional or longitudinal motion 
£ = length of one-dimensional subsystem 


with Eqs 


(41) or (55) and (62), Eg. (63) becomes 


2 


(64) 


' 1,2 


Re(Zj) ReCZ^) 


7Tu»n . 


12 


- mode density of source subsystem 


Eq. (64) IS an a convenient form for evaluating the coupling loss 
factor. Expressions for the mode density of the source subsystem 
available for a variety of subsystem types. The input 

Fo? for many common subsystems. 

For example reference [21] presents expressions for rods (longi- 
tudinal motion) and beams and plates for both force and moment 
inputs • 


For subsystems connected along a line such as occurs between 
transmission coefficient is a function of the 
of the incident energy with respect to the junction. The 
tionship between the coupling loss factor and transmission 

involves an average over the incidence 

[ 1 / J t 


plate 

angle 

rela- 

coef- 

angle 


n 


1,2 


c 

R 

2ujA 

s 


<t(0)cos 0 > 


(65) 


where 


< > denotes the average over angle 

is the length of the line junction between 
the subsystems 

0 is the angle of incidence 


is the area of the source subsystem 

Acoustic spaces are coupled to each other and to bounding panel 
K T relationship between the transmission 

factor is ^ between two acoustic spaces and the coupling loss 


n 


r. A 

_ Y Q w 


4u>V 


( 66 ) 
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where 


w 


the area of the wall connecting the spaces 
volume of the source space 

coefficient which is 
(TL) I? transmission loss 


TL = 


-10 Logjo (i) 


( 67 ) 


applicable'^^fjJr detemfnina^^he Panels is directly 

Eqs. ( 66 ) and ( 67 ) ^ couplxng loss factor according to 

Is%ht%cielizirby\^hrr\dY^^^^^^^ -Paces 

a transmission coefficient Thi<; ^rad' opposed to 

coupling between subsystems w^ic^di m ^|Wal for the 

dimensions that describe the subsystem one in the number of 

sional subsystem while an acoustic panel is a two-diraen- 

sions. The coupling loss space involves three dimen- 

related by ^ ® factor and radiation efficiency are 


pc A 


■1,2 u> M, 


rad 


( 68 ) 


^^^[^teristic impedance of the acoustic 

^p ~ area of radiating panel 

Mp = total mass of radiating panel 

given in^^the^ JStion^on Icou^sUc/ae^^^^^^ efficiency of panels is 
procal relationship for couplYno sources. The reci- 

to determine the coupling loss^ facte/ conveniently used 

and a panel: " actor between an acoustic space 


( 69 ) 
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gg yer input evaluation . - The SEA equations describe the oower 
exchange between subsystems of the SEA model. The inputs to the 

extlrnaPsources'^^^A^T-^^^ subsystems tnrough contacts with 

external sources. A direct measurement of the power inputs under 

actual operating conditions is usually very difficult althnnah 

intensLr »"^a^u/ements offer’^SSe 

pressure levels i" vibration Sr 

conditfnnc^ contact points under actual operating 

conditions and then to relate these measured levels to the nnwer 
inputs through calculated or measured input impedances. 

IS* evaluated 

Second ^ ® mode-by-mode basis. 

^ U ^ using input impedances. The impedance 

ancrdata''^%hl analytically or using measureS^^fmpe^ 

ance data. The following material will show that th#^ 

impedance approaches lead to the same result. 

SiglS‘jiIldS''is”'°‘**'‘’^'"°‘*® “PP'^vech, the equation of motion for a 


u. + n.ui.u. + 
1 111 


“i “i = s fi't) 


( 70 ) 


where u . ( t ) 
viscous ^ damping 
fj^(t) is the 


is the modal response, M is the system mass, n 
loss factor, u>. is the 
modal force given l^y 


0 M * 

resonance frequency. 


is a 
and 


f.(t) = Jdx si. t) • 4i.(x) 


( 71 ) 


where s (x,t) is th^ vector stress acting on the svstem frnm 
external sources, and *l»i(x) is the mode shape, and the intearal is 
over the spatial exter^t of the system. in the case where the 
force IS broadband random stationary time history with a 
flat spectral density, the power input to mode i is qiw^^by 


y(i)_ 2 _[i 
in 2 M 


( 72 ) 


spectral density. The time-average power input 
to a group of modes is found by summing over the group 


( 73 ) 


W. 

in 



S 


f . 
1 


Equation 73 can be rewritten as 


W. 

in 


- 2 N ^ 

■ 2 M f. 

1 


(74) 


where N is the number of inodes in the qroup and ^f. 
over the group. i 


is an average 


If an ensemble of systems is defined in which the resonance 
frequencies are distributed randomly, the number of modes, N, in a 

frei^ency band, Aui, can be written in terms of the modal density, 
n(ut ) , as ' 


N = n(uj) Au) 

The power input in that frequency band is then given by 


(75) 


y(u»)_ n n(u)) 
in 2 M 


(u») Au) 
i 


(76) 


where w is the band center frequency. 


For the case where the excitation is a point force 
of the modal force becomes 


the spectrum 




(77) 


location at which the force is applied and S, is 
becoi^^^^^^™ applied force. The average over the gfoup 




( 78 ) 
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where »l»| is the value of the mode shape squared at location x 
averaged over all modes in the group. 

In systems where the mass density is uniform, the value of IJi? 
tends towards one as the number of modes in the group increases^ 
However, at points near boundaries to the system, the value of ip? 
may be significantly above or below one. ^ 

If the concept of an ensemble of systems is expanded to include 
systems in which the point of application of the force is a random 
variable over the spatial extent of the system, it can be formally 
shown that the average of j|»? (Xq) over the ensemble is one. In 
cases where the mode shapes ^re not known, this result is useful 
and is commonly a part of a SEA model. The power input for this 
case can be written from Eq. (76) as 


w. 

in 


_ n n(u)) _ . 
- 2 M 


(79) 


The formulation above can be expanded to include cases where the 
excitation includes moments and multiple excitation points. It 
can also include a distributed excitation along a line or over a 
surface. For example, the formulation can be expanded to include 
a distributed excitation due to a turbulent boundary layer or 
acoustic pressure field. 

For the case in which a two dimensional structure is excited by a 
distributed pressure field, the spectrum of the modal force 
becomes 


Sf = / dxj 

i 


J dx- S . 

-2 p(Xj) p(.X2) 


4*i(Xi) «P-(X2) 


( 80 ) 


S 

where p(Xt ) p(X 2 > is the cross -spectrum of the pressures at Xx 
and <2 and the integrals are over the spatial extent of the struc- 
ture. If, in addition, the excitation is homogeneous in space, 
the cross-spectrum depends only on the difference Xj - X 2 and 
Equation (80) can be written 

= / dk Sp(w,k) |»i».(k)|^ ( 81 ) 
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is the cross-spectru„, h 

in^fV’" particularly 

average can be taken over the power input or an 

cal e^tf^ force can be used in Equation spectrum of 
timate of the power input. ) to ootain a statisti- 

g^Me"rtie°^“c1i“,^-' de^on- 

frequency, pure-tone force wii-k ^ point force. When a sinale- 
point x^, the response velocity f V,"" de^kt^^ applied 


V = 


F 

2(u»,x ) 
— o 


(82) 


power input Vs Point impedance. 


The time-average 


w. 

in 


= I Re 


(83) 


where Re signifies 
(83) yields 


"the real part of”. 


w. 

in 



|V|2 Re Z(u),x ) 
~0 


Combining Eqs. (82) and 

(84) 


or 


W. 

in 


=7 Ifl" Re 


Z(u»,x ) 
-o 


(85) 


wiabl'es i^^"Eqt“°84rand ^ esl® can'^b"” ' and velocity 

the power input in a band of freguenc!er?i*if g??en‘’lf“ 
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where the integral is over the band Au>. 

The variables in Eg. (86) can often be measured. The easiest mea- 
surement involves the velocity. Measurement of forces requires 
that a force gauge or load cell be inserted between the source and 
subsystem during operation. This may not be feasible as, for 
example, in the case of the gearbox connection to the airframe of 
a helicopter. 

The input impedance of a structure can be determined from measure- 
ments on the structure through use of an impedance head. This is 
a device which attaches to the structure and measures simul- 
taneously both the applied force and acceleration (velocity) at 
the contact point. A shaker is used in conjunction with the 
impedance head to generate the applied force. It is usually not 
necessary for the level of the excitation to equal that of the 
actual source, since the assumption of dynamic linearity is almost 
always valid. However, in some cases it is necessary to apply the 
same static load to the structure so that the measured impedance 
is equal to that seen under operating conditions. 

The voltage excitation to the shaker can be in the form of a 
sinusoid at discrete frequencies or a broadband random noise 
signal. In the former case the sine wave is swept over the 
desired frequency band. The relative amplitudes and phase of the 
forces and acceleration signals are monitored to give the desired 
impedance. For the broadband excitation the outputs are processed 
digitally by sampling the random signals and then evaluating the 
impedance in terms of Fourier transforms of the sampled signals. 

The power flow is dependent on only the real component of the 
impedance or mobility (inverse of impedance). In cases where the 
reactive or imaginary component of impedance is large compared to 
the real component, the accuracy of the measured value of the real 
component is reduced. 

In the absence of accurate impedance measuruments , analytical 
expressions for the input impedance can be used with measured 
velocities (or forces) to determine the input power. 



IT soTv“V"tof Calculated using a modal formulation or 

able^^*1itow'ever?’'*M is^thi**^^°” of the impedance becomes formid- 

^Iculate power input, a statiltical appr^dh^^n 1?^ is 

known from the theory of room acoustics that the real cart of the 

s^«e"iS*a‘fa'riS’^°“' impedance) for a point ^oluL^eflcUy 
source in a large room is equal to the admittance of an infinite 

acoustic space. The relationship between the real oart of the 

eJ?sS"7or "fnv" dv for anTnfiri'te sysS^ 

exists for any dynamic system under the following conditions: 

frequen1?Y, Tr‘*®^ contribute to the response at a single 

® frequencies contain- 

mg inany resonances , or 

IxcftadTom possible locations of the 

With this result the power input from Eq. (86) can be written 


W. = S^(u>) Re 


Auj 


( 88 ) 


aAS^z^^ the impedance of the infinite system, provided that 
b^(uj) and constant over the band Aw. 

Although ^ Eq. (88) is identical to 

prate, the impedance is given by [ 22 ] i-j-ai- 


^inf. plate = 


(89) 


where m^ is the surface density (mass per unit area) k is the 
The Saflenlit/ /o"r“a'^'fTat pllte^i%r>ln'br‘‘“"'"*' “'"*""*""- 


n(uj) = 7-^ 
4hkc, 


(90) 
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a 


where A is the area of the plate. It is left as an exercise for 
the reader to show that Eqs. (79) and (88) are equal with these 
expressions . 

An extension of the impedance approach to cases in which the 
excitation includes moments and/or excitation by coherent sources 
at multiple points becomes a good deal more complicated. However, 
use of the statistical concepts introduced earlier makes these 
extensions possible and many special cases have been worked out. 

Up to this point, the excitation has been a stationary random 
process with a flat or white spectrum. In cases where the spec- 
trum is broadband and varies slowly over frequency, a white 
spectrum gives a good approximation to the time-average power 
input to a single mode as long as the spectrum is flat over the 
resonant response bandwidth. 

If the excitation is at a single frequency, the time-average power 
input can be derived from Equation (32) as 




in 


M 


2 

u) n .u). 
1 1 


.2 2.2 ^ 222 
(u). - uj) +u)n.ut. 
1 11 


(91) 


where F. is the complex amplitude of the modal force and u> is the 
radian rrequeiicy of the excitation. The power input to a group of 
modes is found by summing over the modes . 


w 

in 


= I 

in 

1 


(92) 


Equation (92) can be evaluated exactly in cases where the reson- 
ance frequencies and damping loss factors are known. However, if 
the modes are coupled to modes of other systems, which is usually 
the case, the resonance frequencies are shifted and the damping 
loss factors must be changed to account for the apparent damping 
due to energy transmitted to other modes. Evaluation of the 
resonance frequencies and damping loss factors in the coupled 
system requires solution of the entire system equations of motion 
which negates the simplicity of the SEA model. In addition, the 
accuracy of the computed resonance frequencies and loss factors is 
often poor because the system caniiot be described with sufficient 
detail . 


sys\\^™s‘arrV£^e°d^tn ^wMch^ 'f ^n ensemble of 

domly distributed with a ^^onance frequencies are ran- 

power input to a single »ode i av?ra|eTl^‘i“L 


wfi^= 2 i i ,p ,2 
2 M 2 


in 
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similar Equa^on'^^(^ 72 ^r wiSf i- ensemble. This equation is 
force replaced by the amolitude random modal 

force. This result allow? an sea 

excitation. allows an SEA model to be used for pure tone 

m2n?r?r ti? eLrgLr^f '^he^'^diri measure- 

actual operating conditions ^ excited subsystem under 

power involves ^a meatu^erent of^pplLf f^rce'or' resuf 
velocity, it is oft^^n ;^q • lorce or resulting point 

energy l4vels^!„“t^^rd??ec=??reiiciLrsSi^^Lr^^“^* 

pressure levels*^in*^the ”spac^ surround^na'^th'^'"^^"K *^° the 

measuring vibration levels on gearbox as opposed to 

radiated acous?iJ power is 'torn which the 

radiation efficiency where sir" estimate of the 

involve a single subsvstem i-h#» attachments unambiguously 

system may more appropriatelv response of that sub- 

suring thl attachefel.tTo“nt 4loc^^^^^^ ’^*>an mea- 

power based on the modal or impedan« apprlaches «®timating the 

leiceeie^e" a«?ee“''«speese^''ieeeie^\ subsections the measured 
subsystem energy ® determine the 

give^ byV acoustical subsystems the total energy is 


= V.-D 


( 94 ) 


Where V is the volume of the space and D is the energy density 

P = <P^>, /(pen 

and <p2>jj IS the space average mean square pressure in the space. 
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For structural subsystems 


pStr 

^tot 


M <v^> 

X 


( 96 ) 


av^ge" - «>e space 

eneraies terms proportional from measured 

righr\%„^rs"fde‘^T/ij;i 

solvrng for the system^rlspXf l°evels“'‘ =°«ce"tt™s "i^ 

-lution Of the sea mat ■ 

model. The total enerav fn^« subsystems of 

the modal energy and mode density "‘^Slaing^lo^^"'’^^ obtained from 

F ^ 


E. = ( 
i.tot 


E. 

1 . tot 

n . 


-) 


n. (u») 


(97) 


Where the guantity i„ parenthesis is the modal energy 
The total . 

X/tfms"^"’ s’?/rtu°rl/‘'a^^i-aLS??c1r 

Acoustical subsystems: 

<p 2 > = (futot^ 

IX n ^ 


IL>^ O 

2n^c 


(98) 


wnere <p?) So a-u_ ... 


( 99 ) 


Structural subsystems: 


<v?> = (• 


E. 


1 X 


1 . tot 

n . 

1 


n . 
1 


) 


M. 

1 


density for structural subsystem 

“ total mass of subsystem 

~ space average mean square velocity 

motions involve translational 
plates ^ bending deformation in beams and 


velociuJs‘*®^°™'^"® "" “9dlar 

>x: 


<0,2> - f_ l itOt ' 

IV ^ ^ 


n . 
1 


n. 

1 


p.L.I 

1 1 p,i 


(100) 


section Vlf moment of inertia of the cross 

section, (t. is the density, and is the length of the beam. 

a^Jordfurt^f displacement levels are determined 


<a?> = lu ^ <V. 2> 

IX O IX 


and 


<d . 2> 

1 X 


<V.2> /ui 
IX o 


( 101 ) 


2 


( 102 ) 


Wq = center frequency of band 


” space average mean square acceleration level 
i X space average mean square displacement level 

ments!“' ^or angular accelarations and displace- 
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SEA Model of the Sikorsky S-76 


described the general features and reouired 
steps for implementing a SEA model of noise and vibration trans 
mission. The fundamental approach of SEA was presS in 

the^ve/a°lT'^modil®"'*c*’“®^^°'’® between individual subsystems of 

General procedures for computing the nara- 
meters in the system equations were also described. ^ 


Though the basic approach of SEA is conceptually straiahtforward 
and easily understood its implementation retires a well deveWd 

dynamic behavior of mechanical and acoustical 
systems, in general, and also of the particular structure under 
investigation. Sectioning the structure into multiply-connected 
structural subsections requires an understanding of ?^?irof 

t°o^“Sr"" iunltio„s“iUTo=up?f 


The foi^owing sections illustrate the implementation of the sea 

SrafcirL helicopter that is representative^! Ihe 

g neral ^l^ss of civilian/commercial helicopters. The helicooter 
chosen, the Sikorsky S-76, is a twin turbine powered helicontlJ 
a certitied maximum speed of 80 m/s, a cruise speed of 7S 

wetghro/«",2f “"rrr ^ certfffed gross 

ignr. oi 45 KN. In a standard configuration the S-76 can armin- 

exeJutive comfLt^^^®”^®''^ ^ ® passengers in wide body 


The development of an SEA model for a particular apnlication 
involves a series of steps that are identified in Figure 19 JSe 

desJ?ibina®^SFA <^iscussed at length in the general section 

^ ^ PJ^®sented here in block diagram form to 

help focus the discussion of the SEA model for the S-76. 


The first part of this section involves a description of the s- 7 fi 

SEA^subsvJtemf .^^''siderations leading to the selection of the 
tiin ftf ^ model. A second part deals with the evalua- 

and model including the mode density 

and damping and coupling loss factors. Reference is made to 
Standard expressions from the general section on SEA that are 

r'e^\!/d"rh'^v noh-s?andArd%ressi!S! 

described within this section. An example 
u!ed !x?«!T«iy*Srihe 11-“® honeycomb panels that are 





An important part of this section deals with an evaluatinn of 

tXl “f'SLV'VTanS“*^a°„“e\ 

menbers at junctions with cross frames. A derfvation of ^ 
A^Sil b\°®® eapressions for frame junctions is given in 

Dower^%\ma'ti^?^of evaluation of input 

the sub^-^tem eniroi^t “atux equations, and solution for 

evlluluir Vriie? fK ■ "SP?"®* levels. The input power 

helicopters. The estimation'^ of source*'lwell i s ** dependent^ or i- 

the^ system procedures for formation of 

tne system matrix equations and solution for the subsvstem 

se«?on" discussed in the geneta?®^” 

^gscription of the S-76 ai rframe . - As previously described the 

?l^«Lon *>®4®°P^®-^ ®®l>in noile in J untreated' c^! 

•w qearbox. It is a source of both acoustic anH 

vibratory energy that is transmitted into the cabin acoustically 
through intervening panels and spaces and vibr^ionaTlTthroua^ 

surfaces that raTate dSec?ly 

to thrairSime S-76 is pad mounted directly 

tL aLrb^x attachment points. On other helicopters 

strut? ^ similarly mounted or it may be supported by 

struts and torque restraint itiechanisms . ^ 

Fiau ^?7 airframe structure are displayed in 

overhead framing in the cabin consists of two main 

Difot^Q^w reach from the gearbox forward to the 

pilot s windows. The gearbox is mounted directly onto the fore/ 

inteqJal^V underneath the gearbox the frames are an 

integral I beam construction with added platina for additional 

the gearbox the framing is 2?ivetef ug^^^^ 
>ight construction where the weight decreases moving forward. 
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Identify Subsystems of the SEA Model 


Iden.-fy Types of Motion 
Bending, Torsion, 
Acoustic Modes 


Identify Physical Junctions 
Connecting the Subsystems 
Point, Line, Area 



Evaluate SEA Parameters 


Mode Density 

Plates 

Frames 

Acoustic Spaces 


Damping Loss Factor Coupling Loss Factor 


Added : Damping , 

Absorption 


Frame Junction 
Frame/Panel 
Panel/Acoustic Space 


Power Balance Equations Evaluate Input 

Power 


SEA System Matrix 


Source Levels 


Invert System Matrix 


Modal Energies 


Compute System Response Levels 


Vibration Levels 
Pressure Levels 


Figure 19. Block Diagram for SEA System Response Evaluati. 











and extend down around the sidls th? riveted construction 

structure. The skin panels iorboth ?L^^^\^® support the floor 

® light weight honeycomb construrtf*'®^^ 
the frames. m t-Jie bare condition riveted to 

framing and honeycomb skin panels are treatment the 

t-aneis are completely exposed. 

4-U^ _ J >1 


between the^ overheld^^fraSies^^ ^ frame structure 
the broom closet**) frames and the floor (call^^fi 

tail rotor controli ' pa 2 "‘, 5 rrouof the JL“nd 

overhead surface of the cab °n t?' s structure to the exterior 
helicopter extending forward from section of the 

pilot/copilot 's seats is a fih<»rrr? -raming section over the 
section is a typical rib%tr\ngi^%"?^|=cVsrrS°"- "*'* 

'TKiOt 


-• -r - W..1. Wiw A WIl * 

Which ^ alio ^covers^5ie^®co^^ ® fiberglass cowling 

skin surfaces forward of rh<» inkages mounted to the exterior 

cha^Ir’Sfthe turb?i!?'*^1 separate's^tee 

The plenum chamber is also seoarat-oH f inlet air supply, 

by a bulkhead. Immediately aft of the^rea^^^ turbines themselves 
a luggage compartment that lies blllw rho is 

inlet spaces. Behind it and A*'® gearbox and turbine 

occupied by the envi?oii,ental contro^^ ^^^^i«®s is a spact 
cabin are shown 'in Figure 2 ” control units. Dimensions of the 


in a 

tions of the total structure which Ixhi ®®c- 

dynamic behavior is chlJ^iri^- dynamic be- 

deformation within the structurl af^t /11 type of 

coupling to adjacent structurac d nature of the 

include bending or flexural tnrci types of deformation 

motions. Each differs in thV wa? in° 2 hJ;>."°"®^*“‘*‘'*®^' anS shear 

the structure. ^ which energy is stored within 


b"^®d coup^^^ motions 

d!??lrenf V°Vs‘' o?"'o’^*l®' n®vV co,iS !Sf“"««®"thri.rd*eri!2b 
bSd^treated collectively as a S?S|i®|®"s'4sTstei ^ 


section the^*total ^structure ^nto^*ind'^® physically 
Structural subsections. This ie J^ndividual components called 
considerin, the physical ^ 


(7r 
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ORJG(r;/?L PAV.V.: XI} 
OF POOR QUAulTV 


sectioned by the framing to which they 
fJLes This is h® sectioned by junctions with cross 

llllr 1ife° 

def^rSise'^space?.' surfaces which physically 

In many situations on the actual structure there are potentiallv 
structural non-uniformities in the form of cutouts^ 
cabling, hydraulic lines, etc., stiffening 
a2oiS?T-v ^ plates, as well as local changes in cross-sectional 
2i?h ff acoustic spaces the partitions may be only partial 

Witt significant open areas between the spaces. Decisions must be 
made concerning the degree of detail that is accounted for in 

subsections. Accounting for additional detail 
results in a larger number of subsystems. The costs to implement 
model increase, and there is also the potential for 
greater computational error. The trade off must also take into 

accuracy of the modeling procedure itself 
...n deciding the degree of detail included in the model. 

Ivl subsectioning of the S-76 airframe for the initial 

between^ fri^es^^'^®frA*l ^ includes the panel sections 

window ' ^ ^ sections between junctions, bulkheads, 

window and door panels, and assorted otter structural elements 

cabin itself other acoustic spaces include the 
different overhead compartments around the gearbox, turbines the 

iisted^in Tabir‘ 2 .^”'^ luggage and ECU compartments. Thes4 are 

®®^^ structural subsection, subsystems have been included 
motion groups Containing different types of 

o?^o K primarily for the framing at lower frequen- 

where tte final subsystem selection accounts for torsional 
motions and bending motions about two axes perpendicular to the 

nJlt® model as formulated for the S-76 does 

longitudinal motions in the frame due to the 

attattmJit motions, which predominantly excite the 

attacnment frames in bending and torsion. 

c 7 A Structural subsections are used in describing the 

S-76 airframe . Of these 45 are frame structures and 37 are panel 

ySwinrV ‘hree types of motion are coS!dl«d 

Ik total of 135 fr^e subsystems in the SEA model. In 
addition there are 8 acoustic spaces listed in Table 2 that are 
included as subsystems in the SEA model. This brings tte total 
number of SEA subsystems for the s-76 to 180. tte total 


0 
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environment, The^^oSfe? spaces interior noise 
turbine spaces) or as co£?ectina (e- 9 -/ gearbox and 

ment). Others, such as th#» nr>^A ^ t®*g** luggage compart- 

control unit space, are less environmental 

they physicall? are not on a Xecf ^Paces since 

cabin. They are of importance ii ^ i" source to the 

distribution Of aooustL levSfs^tt^ugho^ftKrKpLr"""" 

the airframe about^ h^^for^aft *center*line'^**Th'^** symmetry in 
reflected in the structural siAseot^f^nt ®’ synroetry is also 

or frame sections on Ue left In^r^aht 

identical surrounding structural have 

for the different lourSf are af/n i«P^ts 

energies of the symmetric subsystems wiir“h^''Jh' 
allows for potentially sionif^ranr P®. same. This fact 

and solving the power Llance equaUons?^^^'^^^^®''® formulating 

ant m^oTle„^\ir"are*rd^^^^ 

degrees of freedoT (il ilZlJ the number of 
system matrix Th^ effict on energies) in the SEA 

no symmetric counterparts is \o corS^/nif ®^®y®tems with 

terms for the symmetr^ pairs off-diagonal coupling 

s?rucSrl s";Te\"rx/ "^® 

smaller than the original total of significantly 

subsystems accounting for bendina total number of 

frames is 106 as comp^arf^to Tso oritnally 

depend not onlt* n the Sure^ofth these reductions 

of the excitation arwellT ?? 5 , ^ also 

to excitation at only one of fhl considering the response due 

the problem is not symmetric Also Vf points, then 

tions are not left/riaht actual gearbox excita- 

symmetric and must be solved ^ its full 

high deT«V\f®VnTeVco1mlctedness®h'’V^^ airframe there is a 
sources %o receivers ?f*s diflw^^ P^^hs from 

transmission in terms of nathc i-Kn-h describe the overall 

in terms of paths that are independent of each other. 


93 




The extent to which this is possible depends on the existence of a 
se<^uence of strongly coupled subsystems with relatively larger 
modal energies than surrounding subsystems. in this case the 
power flow will be along the path of the strongly coupled sub- 
sy.^tems with lateral flow out into the surrounding subsystems. 

Prior measurements of a preliminary nature on the S-76 have 
suggested that the flow of vibratory energy from the gearbox along 
the main fore/aft frame members constitutes an important path in 
important frequency bands below 1 kHz. The gearbox is directly 
and strongly coupled to this path and as well it is anticipated 
that the coupling between sections along che framing is strong. 
Important lateral flows of energy from the path would be into the 
adjacent skin panel sections and also the cross frames at frame 
junctions . 

Other potentially important paths from the gearbox that are 
accounted for in the SEA model are into the drip pan immediately 
beneath the gearbox and the rear passenger bulkhead. Both are 
excited by vibration transmission through the frame to which the 
gearbox is directly mounted. 

Acoustic spaces in the helicopter airframe are an integral part of 
the overall SEA model. Energy stored in the resonant acoustic 
modes of the space are coupled to the resonant modes of adjacent 
structural subsystems. In addition the spaces are coupled to each 
®^her directly through non-resonant mass-law controlled response 
hhe intervening panel structures . The presence of holes or 
openings in the panels also contributes to the direct coup- 
ling between acoustic spaces. 

The number of acoustic spaces in the helicopter airframe is small 
compared to the number of structural subsystems. In some cases 
the transmission is from the source space through a panel struc- 
hure directly into the cabin such as occurs for the gearbox 
acoustic space through the overhead panels into the cabin. For 
other paths the transmission may involve a single intermediate 
acoustic space. Airborne transmission from the spaces around both 
the gearbox and turbine can involve the luggage compartment as an 
intermediate step in the path to the cabin. 

Vibration and acoustic transmission in the helicopter airframe are 
inherently coupled. The degree of the coupling depends on the 
parameters of the particular model so that it may be appropriate 
to describe independent airborne and structureborne transmission 
problems. For helicopter airframes both probably must be con- 
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Table 1 - Structural Subsections for S-76 SEA Model (cont'd) 
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honeycomb panels and fiberglass 


Table 2 


Acoustical Subsections for the S-76 SEA Model 
D escription Identification Code 


Space around gearbox and under forward 
cowling 

GBA 

Turbine inlet air plenum 

TIA 

Space around turbines, left and right sides 

TAL,TAR 

Luggage compartment 

LCA 

Environmental control unit space 

ECUA 

Cabin 

CA 

Nose compartment 

NCA 


99 




Mode density for the S-76. - The S-76 model is comorised of fr;,«>*. 
iSin. ’^*’® applicable equations are shovm in 
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Table 3 


Summary of Parameter Requirements for Frame and 
Acoustical Subsystems 


Frame Subsystems 
Torsional resonances 

Bending resonances 

S-76 frame substructure 
#'s from Table 1 


Applicabl e Equations and S-76 Subsections 

Eq. (41) with wavespeed from 
Eq. (43) 

Eq. (53) with wavespeed from Eqs. 

(49) and (50) which yields Eq. 

(55) ^ 

1 - 32 , 51 , 52 , 61 , 62 , 70-78 


Input Parameters Required 
- length of frame subsystem 
Torsional Motion 


Symbol 

i 


- torsional moment of rigidity j 

- polar moment of inertia of I 

cross-section P 

frame^axisl^°” (about each of two axes perpendicular to the 


- radius of gyration of 
cross-section 

Material Properties 

- Young's modulus 

- Shear modulus 

- density 

Acoustical Subsystems 

Acoustic resonances 

5—76 acoustical subsystems 
from Table 2 


K 1 , K 2 


E 

G 

P 

Eq. (40) 
Total # = 8 
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Table 3 (continued) 


Input Parameters Required 

- Volume of space 
“ Surface area 

- Edge length 
Physi cal Properties 


- Speed of sound - (Possible 
temperature and ambient 
pressure effects) 





V 

S 

P 



For plate sxibsys terns additional expressions are required charac- 

Conveitfon^?^ dyn^ic behavior of composite honeycomb panels. 
Conventional single layer panel subsystems include the rear 

and bulkhead, windows, the drip pan underneath the gearbox, 

surfaces on the doors and nose cone. For the conven- 

aia?vtin layer panels the required parameters and applicable 

analytic expressions are given in Table 4, 
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Table 4 


Summary of Parameter Requirements for Single 
Layer Panel Subsystems 


Panel Subsystems 


Applicable Equations and S-76 
Subsections 


Bending resonances 


S-76 Single Layer Panel 
Subsystems from Table 1 


Egs. (49), (50), and (54) which 
yield Eg. (56) 

35, 44, 45, 46, 49, 50, 53-56, 
59, 60, 63-66, 69, 79-82 


Input Parameters Recmired sy«>h«^1 

- Area of panel A 

Thickness of panel h (used to compute 

K , radius of 
gyration) 

Material Properties 

- Young's modulus 

- Poisson's ratio 

- density 


0 



/ 


» iUustrated 

three distinct regions. involves 
the honeycomb core material i?^ sumcie^^^^^ panels 
tively prohibit any thickness deforma?i^<^^^^? ef fee- 
very high frequencies, well above th?^range of 

tional to the^^^^are^root wavespeed increases propor- 
for a panel def^Sn^in the reLlt 
composite cross-section i <5 i- tending rigidity of the 
in the face sheets Vhi\ art extensional deformation 
honeycomb core. paced apart by the thickness of the 


2(l-a2) 


( 103 ) 


The expression for the wa\'espeed is : 


thickness of face sheet 

- thickness of honeycomb 
core 

= Young’s modulus of face 
sheet material 

= Poisson's ratio of face 
sheet material 


( 104 ) 


Pfs “ surface density of one 
face sheet 

Pc ~ density of honeycomb 

core (in expanded form) 

bute signifiSa?Uy''to''throvera^^^^ tigidi?^ . 

important /7nd‘*th^ *«aves|«r Vn become 

independent of frequency consistenHi?!; shearr.tl pTopa,SI?oi'' 


0 



Figure 21. Dynamic Behavior of Composite Honeycomb Panel 



The value of the vevespeed for sheer wave propagation is: 




( 105 ) 


G - shear stiffness of 
honeycomb core (in 
expanded form ) 


deformatio^n^bSt wTth^^a'^troalle? efflct°^ represents bending 

bending rigidity in Regions Ts dlt^rli?^^^ ThI 

sheet deforming in bending: determined by an individual face 


Et| 

D, = £2— 

3 l2(i-o2) 


(106) 


The wavespeed is given by: 


(Pf,*iP,H) 

D 


(107) 


frequency fro^Egs^^UMfraoir o/TlVf and® 

In on 1 • ^ 


In Region 1 : Low frequency bending of the entire composite 


rij(w) 


4n (tt 1 1- 

p H)' 


(108) 




In Region 2; Shear deformation of the honeycomb core 


02 (w) = 


Am 
2n c| 


(109) 


In Region 3: High frequency bending of the face sheets 


113 (m) 


An I 


D. 


(P 


fs 


+ 1 pH) 
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S“-f 

gent section, the radietiofimlienc? “aj be tlh« tS be^^^ 
within "the yanel^.* <^>>e type of deformation occurring 

si a^soL^ty tlV p° ^f^W|:76 a-^o/" structures 

raJgrof vlluU"ney"' structures d^ping^hatt Tbroad 

manufacturers' estimates* o^f the damping materials 

values will be used resulting composite loss factor 
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I I aluminum honeycomb 
tiv-H riBERCLASS honeycomb 


Figure 22 S-76 Plate and Skin Materials 






panels are either of a honeycomb or fiberglass construction. 
Damping values for the fiberglass material itself depend on the 
formulation, the type of resin used, the density and orientation 
° Reported values are in the range from .003 to 

Ki Additional information will need to be obtained, if avail- 
able from the manufacturers of the fiberglass sections. A 
question for the built up structure is whether there are signi- 
f leant contributions from interface damping due to gas pumping or 

itself The greater the inherent damping of the material 

Itself the less will be the influence of the other sources of 


panels the damping levels will exceed the 
levels for the individual component materials alone. The com- 

inciudiia^tht depend on the nature and quality of the bonding 

in5irc adhesive used. In the absence of manufac- 

K ^ described in a subsequent section, 

needed to estimate damping levels for the honeycomb 
panels. Particular front window designs involve a multi-layer 
construction with an adhesive middle layer. Here damping levels 

estimated based on experimental measurements on the 
acuuai niacenal* 


A difficulty in performing such measurements is isolating the 

energy which is transmitted to 
interface damping is important. 
It IS not feasible to simply remove the panel of interest from the 
remaining structure since this also removes the desired contribu- 
to interface damping. For structures where interface 
damping is not expected to make a significant contribution as for 

panel constructions , the panel may be tested separate 
rrom the attached structure . ^ 


Damping for the acoustic spaces in the S-76 model includes only 
dissipated within the space and not the energy trans- 

viS^atinn° coupled into resonant 

yibration of the bounding panel surfaces of the space. Air 

small spaces on the helicopter in the frequency 
interest is negligible. Where added treatments provide 
absorption in comparison with the totally bare con- 
figuration the damping can be estimated based on the supplied 

characteristics of the treatment while completely 
neglecting other contributions . ^ ^ 





Estimating damping levels of acoustic spaces based on experimental 
2^ reverberation time is feasible to the extent that 
the dissipated energy within the space is large compared to the 

depends in part on the transmission loss 
panels. Where holes or gaps contribute 
significantly to the transmission, covering them can reduce their 

contribution and improve the validity of the damping estimate for 
uxic spac6 . 

Damping level estimates are a source of uncertainty for all 
vibration prediction methods involving the resonant response of 

1 finite element methods, the methods 

mechanics, and of SEA as well. Damping does not 
tl convenient analytical prediction procedures for 

of Unfortunately it is important for the prediction 

of vibration response and transmission. ^ 

Coaling los_s factors for the S-76 . - The required coupling loss 

model in Tables 1 and 2 are identified by 
the junctions between the subsystems. The junction types for the 
S-76 model include point, line, and area connections. The frame 

subsptems are considered to be connected to each other at a 

Junctions between frames and panels occur along the Ic. gth 
over their^aVeas”^ panels are "connected" to the acoustic spaces 

To illustrate the identification of coupling loss factors consider 
the frame junction shown in Figure 23. The junction is located at 
cross frame location 3 between the rear passenger door and mid 

of cabin. 

Four frame substructures are connected at the junction. Their 

numbers are as follows: ^neir 

Substructure # 

(from Table 1) 


1) 

3 - 

LF23L 

Main longitudinal frame 
frame locations 2 and 3 

2) 

5 - 

LF34L 

Main longitudinal frame 
frame locations 3 and 4 

3) 

26 

- CF3L 

Cross frame at location 
section 

4) 

27 

- CF3M 

Cross frame at location 
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The allowed motions for frame subsystems include bending about two 
axes and torsional motion. In total 12 different SEA subsystems 
^®r^®cted at this one frame junction, with complete genera- 
lity where all of the subsystems are connected with each other 
there are 66 different coupling loss factors characterizing the 
junction. 

In actuality the total number is not this large. As a result of 
the geometry of the junction certain motions are simply not 
coupled to each other. For example, bending motion in the plane 
of the framing is not coupled to either torsional or out-of-plane 
bending motions for any of the four frames . Symmetries occur when 
opposing frame substructures have identical cross-sectional 
geometries which further reduce the total number of coupling loss 
factors. For a symmetric frame junction the total number of 
independent coupling loss factors that must be evaluated is 22. 

further, all of the frames are identical the number of numeri- 
cally different and independent coupling loss factors is reduced 
to 11. 

Frame junctions are the most complicated as a result of the large 
number of subsystems involved. Other junction types are consider- 
ably less complicated. Consider the line junction between a frame 
and the panels on either side to which it is attached. The total 
number of subsystems is five; three for the frame and one for each 
panel. For the completely general case the total munber of 
coupling loss factors is 10. For identical panels on both sides 
of the fr^e the total number of coupling loss factors to be 
evaluated is three. Area junctions between a panel and adjacent 
acoustic spaces are the least complicated and involve only two 
coupling loss factors between the three SEA subsystems. 

Frame junction coupling loss factor evaluation for the S-76 . - The 
analytical procedure for evaluating the coupling loss factors at a 
frame junction is described in Appendix B. The approach involves 
a computation of ener^ transmission coefficients where the frames 
are extended to infinity away from the junction. Reverberation in 
the frame sections is accounted for in the expression relating the 
transmission coefficient to the coupling loss factor. 

Frame-p anel coupling loss factor evaluation for the S-76 . - Frame 
subsystems are connected to panels along the length of the frame. 
The subsystems di ffer by one in the number of dimensions that 
characterize them. The frames are one dimensional along their 
length while the panels are two dimensional over an area. This 
situation is more appropriately described by a radiation effi- 
ciency for the one dimensional subsystem, i.e., the frcune, which 
relates the vibratory motion of the frame to the power it radiates 
into the plate. 
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Analogous to the case of radiation by a panel into an acoustic 
space the radiation depends on the wavespeed in the frame relative 
to the panel wavespeed. When the frame wavespeed exceeds the 
value for the panel a matching or coincidence condition occurs 
where power is radiated more efficiently into the panel. The 
behavior is more complicated than for panel and acoustic spaces 
because even above coincidence a near field disturbance remains in 
the panel due to the nature of bending deformation. 

Below coincidence for an infinitely long connection no power is 
radiated into the panel, and the radiation efficiency is zero. 
For finite length connections power radiates into the panel as a 
result of the discontinuities at the ends of the frame. Un- 

fortunately expressions have not been derived to predict the 
radiation efficiency for the different frame motions below coin- 
cidence. 

The approach followed here is to evaluate the radiated power for 
uniform motion of the frame along its length. This corresponds to 
normal radiation at right angles to the line junction. For 

torsional motion of the frame, 6,, the total radiated power is 
given by the following: 


■1 H = 

rad 

where Z 


L {«ilRe(ZJa2 + Re(Z^)| + Re(Z^ : ] } |6^|2 

I m ni , rj I 

^ = transverse force impedance of panel 

= moment impedance of panel 


( 111 ) 


Z^ ■ = cross coupling impedance between moments 
''' transverse displacements forces angular 

rotations 

a = distance between axis of torsional motion 

and line connection to panel (half width of 
frame where ir connects to panel ) 

L = total length of line junction 

Expressions for the impedances are contained in Table 7 of Ap- 
pendi.x B. In the equations of Table 7, B is the bending rigidity 
of the panel per unit length along the line connection. For 
composite honeycomb panels at low frequency the expression for Dj 
in Eq. (103) is substituted for the bending rigidity B. 
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The radiated power is related to the coupling loss factor and 
frame motion according to the definition of the coupling loss 
factor . 


t ft 

rad f,p tot,t 


( 112 ) 


tot, f 


total energy of frame 
in torsion 

frame to panel coupling 
loss factor for tor- 
sional motion in the 
frame 


where 


'Lot , f 


L p . I , 
p,f 


( 113 ) 


I = polar m<^ment of inertia 
P' ^ of frame cross section 

= density of frame 

These three equations, when combined, give the coupling 
factor between torsional motion in the frame and bending motion i 
the panel as follows: 


'f,p ujp 


(| (Re(Zj)a‘ + Re(Z^)l + Re(Z^ ,,,)a) 

f p,f 


(llA) 


For bending motion of the frame, Oj* the following expression for 
the radiated power results: 


^rad ' - 


( 115 ) 


The radiated power, coupling loss factor, and amplitude of the 
transverse bending motion are related by: 
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( 116 ) 


n 


rad 



E 


b 

tot,f 


'tot, f ~ frame energy in 

bending 


and 


'f,P 


coupling loss factor 
between a frame and panel 
for out-of-plane bend- 
ing of the frame 


llfl' 

~ total mass of frame 

The coupling loss factor is given by 


( 117 ) 


f.P 2wM, 


Re(ZJ 


( 118 ) 


f^Le^'Td^t^^chld'paL ^--tors between a 

tie approximations male In vina'^thell"^’^®' of 

assessed in Phase li by desianinJ f expressions will be 

tie coupiin. loss'' ^s^lol^l^^p^ 

^ -acoustic space no„p,in^ ,n..^^P evaluation tor eK. . .. 

space has been described in generar“'lt^?l ha" adjacent acoustic 
radiation efficiency for the ninii' u 5 evaluating a 

relative to the coilicflince \re|“icy I" ^teguelcy 

the determination of the coincidenr^i single layer panels 

efficiency is straightforward r see^ ^^®^wency and radiation 
dynamic sources for expressions for the acoustic/aero- 

radiation efficiency) As Hicr-noe a ^^^^cxdence frequency and 
composite honeyco^Une’s d/sn^ this section, 

regions. This potentially affeSi the^/^ ®J^^ behavior in three 
efficiency on frequency! ^ tfects the dependence of the radiation 
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1 -hA^ behavior involves bending deformation where 

the bending wavespeed is given by Eq. (104). The coincidence 
frequency for this type of deformation is determined by settina c 

solving for u. (i.e. 2nf T R^dia’ 

are^thin coincidence for deformation ifi Region 1 

re then determined from the above value for f and Eq (6) which 
IS repeated here for convenience ^ <3 1 ) nich 


PA 


rad An^ 


(f)'* 

c 


p 


( 6 ) 


P 

A 


- perimeteor of panel 


P = 


wavelength at coincidence frequency 
(=Vfc) 

boundary condition factor 


P “ 1 simple edge support 
P ~ 1 clamped edge support 

frequency region shear deformation in the core is 
impoiLant, aixu Uie wavespeea is xndependeaL of frequency. If the 
shear wavespeed, as given in Eq. (105), is greater than c then 
conditions are above coincidence and the radiation efficiency is 
set equaJ to one. When the shear wavespeed is less than c , the 
ciency^^^ approach is adopted for computing the radiation^effi- 


The 

can 


radiation efficiency is given as a function of (f/f )^. 
be recast in an equivalent form for bending defor^Tation 


This 
as : 



(119) 


Aq = acoustic wavelength at 
frequency, f 

A = panel wavelength at 
^ frequency, f 
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The 


expression for the radiation 





*^^^iciency now becomes: 


( 121 ) 



speed''of°Eq^^ a05)^“*It°fhoihd‘ b from the shear wave 

Cfency r„ Ee,i'„ 2 ^ari^s 1 re?ae^ “ftl? ®«' 

C 3 , rs fee 

coupUng’\"oVs “rad - then''l«‘7 ^fi, 

couDlinf^ -^aces for the <;-7 a _ ni red 

One involves the acoustic transmission^ sources, 

intervening panel between the spLes openings in th« 

lon-resonant response of the panel ' other is due tc 

accounted for in the SEA model motion is 

adjacent sparL factors betweer 

trolled panel motion that is stronni or mass con- 

acoustic resonances can constitni-o spatially to the 

coupling at frequencies below ^S^cid\'lice"!^°''^^''^ 

acoustic impedance'^Tf ^the^open^^^^ opening depends on tne 

transmission coefficient 


= 


Rp(2^)pc 


( 122 ) 


118 


I 




0 


where 


For circular holes [24 J 




pc = 

\ = 

A = 
P 

2 .. = 


i 8k.j 
iTI 


acoustic impedance 
area of opening 

area of panel between the spaces 

acoustic impedance of opening 
looking into one of the spaces 


(123) 


acoustic wavenumber = ~— 

2nf 


a - radius of hole 


For rectangular holes [24J: 


Z. = pc 

h 16 i)n 


^a£+a_b+^ 
^ a+b 


(12A) 


- -limeisions of hole 


lenqto iJ when t 

arka . ^ compared to the dimension of t 

as ka . o). For circular holes Eq. (122) Secome* 


acoustic wave- 
opening (i.e.: 


= 


.7 


It 

2.88 


(125) 


at low frequencies which reaches a constant value 


of 
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In the other limit when the opening is large compared to the 
acoustic wavelength (i.e.: as ka ^ ») the transmission coeffi- 

the value = A, /A For multiple openings 
the transmission coefficient values” fcPr each opening computed 
sepaptely are simply added together. The calculations are based 
on the same panel area between the spaces. 

Acoustic transmission between the spaces through non-resonant mass 
law controlled motion of the panel is evaluated from the followinq 
expression: ^ 


u>p 

'me = 


-)^J 


( 126 ) 


u) = 2/xf, frequency 


Pp = surface density of the panel 


The result applies for diffuse acoustic fields in the spaces. 

The total direct transmission between the spaces is the sum of 
individual contributions due to openings and non-resonant trans- 

mi 1 on : 


'.1 = ^ 'h,i * 'mJ 
1 


( 127 ) 


where the sum is over all of the holes connecting between the 

spaces. The coupling loss factor is evaluated from and Eq. 
( 66 ) : d ^ 


c A 

n = 

'12 4uiV| 


( 128 ) 


- volume of source space 



Input power evaluation for the S-76 model . - Important sources for 
the S-76 include the mechanical connections between the airframe 
and the hydraulics system and the gearbox as vibratory sources 
where they attach to the airframe and also as an acoustic source 
due to radiation from the gearbox casing and turbulent boundary 
layer pressure fluctuations over the exterior skin. The gearbox 
is the dominant source for both the untreated and treated cabin 
interior . 

The procedures for quantifying the input power for each of the 
above sources are described in previous sections dealing with 
acoustic/aerodynamic and vibratory sources or in the general SEA 
section dealing with the input power. For the TBL pressure 
fluctuations the results in Equations 9-13 and 79 and Figure 17 
are applicable for evaluating the source levels or input power. 

For vibratory sources involving point attachments to the airframe 
Equation 8b is used to evaluate the input power in terms of a 
measured attachment point velocity. Experimental measurements or 
statistical estimates are used for the attachment point impe- 
dances. The measurements are described in a subsequent section 
and the statistical impedance description is discussed in relation 
to Equations 79 and 88. 

Acoustic sources are characterized in Equation 4 by measurements 
of vibration levels, for example on the casing surfaces of the 
gearbox, and inferred values for a radiation efficiency. An 
alternate approach is to measure the acoustic pressure levels in 
adjacent spaces into which the source radiates. This approach is 
described in relation to Equations 94-96. 

The SEA matrix equation for the S-76 model . - The form of the SEA 
matrix equations in terms of the subsystem modal energies, the SEA 
parameters, and the input power is described in a prior section. 
Cambridge Collaborative has developed a proprietary computer code, 
called SEAM, that assembles the matrix equations in terms of the 
SEA parameters and performs the matrix inversion to solve for the 
subsystem modal energies. The program also evaluates the response 
levels according to the procedures described in the general 
section entitled system response evaluation. 

The program evaluates mode densities and coupling loss factors for 
a variety of relatively standard subsystem and junction types. 
Its usage in Phase II to generate predictions for comparison with 
measured data from the flight measurements will involve adapta- 
tions to account for cases particular to the S-76 that have been 
described in this section. 


MEASUREMENT PROGRAM FOR MODEL VERIFICATION 

Objectives 

TaLratory expar!me«rand fn ^ ground/ 

Objective Of and in-f light measurements with the general 

an Phase I. Specitac^ o^hjejtxv'et «e """"foU 

• to provide intermediate levels along transmission nAth<= 

foi comparison with the model predictions ^ 

• to provide cabin interior noise levels and vihrAi--ior> 
levels on panel surfaces that radiate directly JnJf^he 

to provide data necessary to define or validate valuer 
of important parameters in the model (e.g structural 
factor?' ®<^oustic absorption, vibration coupling loss 

to provide data relevant to assessina valiHii-v 
general predictions of the model (e.o. =^„rce ioiere,,?5 

the^®airf?am“jJiS.T’' '^he relative importance' in 

frtmarg'^rihrpaner/i transmissron in the 

ISen?!ff‘:eak“?serTn “"^er experrmentation are intended to 

values /et^^^^r^L^VsL ‘f^o^r rTf^nVng 

SryTes\“h Ve'Vn rtjoht"''” '^""'^oraes: 1 , ground or labora- 

cal shakers or acousUc speakerT‘'fnd® 2 rn’^^^^ mechani- 

with the actual in-fi-in^i- o ' flight test measurements 

dynamic noise sourr#»Q sources (gearbox, hydraulics, aero- 

"bsertions' 1h^“?h“se'"n S.^and otL/"frt^ra:e 

the total program rs sho™ !^^°tji!i 


0 
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The intended end result of the Phase II measurements and com- 
parisons is a cabin noise model that: 

• has been thoroughly validated in terms of the 
fuselage dynamics and acoustic characteristics, and 

• is readily applicable for evaluating the effective- 
ness of a variety of noise control treatments 
(Phase III). 

Ground Test Meaisurement Program 

The ground testing is designed to provide information for validat- 
ing individual sections of the model . The primary focus is for 
vibration transmission in the airframe (see Figure 2), from the 
gearbox attachment locations through the frames and skin panels to 
the panel surfaces that radiate into the cabin. Acoustic trans- 
mission between the different compartmented spaces (i.e. around 
the gearbox, the turbine inlet, luggage compartment, etc.) and the 
cabin will also be investigated. A mechanical shaker or acoustic 
speaker is used as the source. Intermediate levels are measured 
along the path from the source to the cabin where the measured 
levels may be either the acoustic levels in the cabin or vibration 
levels on panels that radiate into the cabin. 

Such experiments are necessary to validate an analytical model fo^ 
a system as complicated as a helicopter airframe. The modeling 
approach. Statistical Energy Analysis (SEA), describes the system 
in terms of the energies of individual sections or subsystems. 
For structural subsystems such as framing members the energ|ies are 
related to average vibration levels on the frames, while for 
acoustic subsystems the energies are related to the average 
pressure levels within the spaces. These extensive measurements 
involve a direct determination of the distribution of subsystem 
energy levels throughout the helicopter for individual source 
locations. Between three and five measurement locations 
required for each subsystem. The SEA model accounts for indivi- 
dual sources for each subsystem and is directly applicable for 
modeling the distribution of subsystem energies in response to an 
individual source. 

Other source locations in addition to the gearbox attachment 
points include locations where the hydraulics lines are attached 
to the airframe. 
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^irfSir„^^H^/rtHe“^\^„^-„i^ rio-/-;;- /„-^y 

uration is intended to sub<?^!:»ni-i -*- latter config- 

the SEA airframe model that preliminary prediction of 

vibration traSission Pa?^ SL constitutes the important 

reader is referred t^TaM? i ? ® gearbox to the cabin. The 

subsections for the Sikorsky ^76^ ^ ^^"‘Pilation of the structural 

gearbox® attactoent pofnts For‘'%hai!f'^®"®® differant 
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Attach ment point impedances . - Attachment point impedances are 
required to appropriately determine vibration source levels 
independently of the structural attachment to the airframe. 
Source levels in the form of free velocities are determined from 
measured in-flight vibration levels and the impedances looking 
into the source and airframe at the attachment location. Although 
the primary interest is with the gearbox attachment points, 
measurements also will be performed for the attachments between 
the hydraulic lines and the airframe. Of particular interest 
be the hydraulics attachment points at the cabin rear bulkhead. 
The impedances will be determined using conventional impedance 
heads and mechanical shakers. 

The shaker is connected to the structure through an impedance 
head. The impedance head provides electrical outputs proportional 
to the force applied to the structure and the acceleration re- 
sponse at the attachment location. The electrical excitation can 
either be a swept sinusoid or a continuous broadband random noise 

signal . 

The impedance is a complex quantity with real and imaginary 
components. The phase angle between the force and acceleration 
is, therefore, important in evaluating the impedance. Also, 
impedance relates force to velocity, and this requires integrating 
the acceleration signal to obtain the velocity. 

For a swept sine wave excitation, processing of the output signals 
is most commonly done on an analog basis, including the integra- 
tion to obtain velocity. If a broad frequency range is in- 

terest this process can be relat.ively slow. With the advent of 
two channel signal processing analyzers, the broadband rundom 
excitation is to be preferred. All of the frequency information 
is generated at the same time. Such instruments are quite flex- 
ible in presenting the results in the form of plots of the real 
and imaginary components or of magnitude and phase. They are 
capable of performing the required integration of the acceleration 
signal as well . 

Dampin g loss factors . - Damping loss factors are extremely im- 
portant in estimating the flow of resonant vibratory energy in a 
structure. For aerospace structures with riveted joints, honey- 
comb panels, and a large amount of attached structure the damping 
levels significantly exceed levels for the base materials. 
Damping levels for built up aerospace structures are typically 
estimated based on general rules of thumb which have been de- 
veloped from past experience. 



126 


For structures where the connections between attached components 
significantly contribute to the damping (so called interface 
damping), it is not feasible to use measured damping levels for 
the unattached components. Measurements on the attached structure 
in addition to actual damping losses include losses associated 
with transmission across the interfaces into adjoining components. 
The total measured loss factor includes damping and coupling loss 
factors . 

Values of total loss factor for different structural subsystems 
determined from decay time measurements on the built up structure 
will be compared with model predictions in assessing the overall 
validity of the damping and coupling loss factor estimates. 

For some structural components it is feasible to measure damping 
loss factors separate from the rest of the structure. This is the 
case where the damping levels of the component are expected to be 
greater than interface damping as a result of applied damping 
treatments. In addition to skin and framing panels with applied 
damping treatments, measurements will be performed on multi-layer 
window configurations with lossy constrained adhesive layers. 

The measurement of damping loss factors involves a straightforward 
experiment in which the steady noise excitation to the structure 
is abruptly turned off and the decay of reverberant energy is 
measured. The energy decay occurs according to the following 
relatioiiship : 


-nuiot 

E(t) = E e u> - center frequency of (129) 

O O . ., .,1 

excitation noise bandwidth 


The measured reverberation time, T , is the time required for the 
energy level to decrease by 60 dE from the initial level. The 
reverberation time is related to the damping (n) by the relation: 

o r 


The energy level of a structure is proportional to the mean square 
vibration level over the surface. It is sufficient to measure and 
average the decay at a limited number of representative locations 
using an accelerometer. Three to five locations are quite ade- 
quate for most situations. 
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Couplinq loss factors. - The determination of appropriate coupling 
loss factor values involves laboratory measurements on component 
sections of the overall structure by focusing on individual 
structural connections or junctions. The experiments result in 
empirical estimations of the coupling loss factor between the 
attached subsystems. One of the structural subsystems is excited 
with a shaker, and the average vibratory energy levels are mea- 
sured in the source and receiving structures . Vibration trans- 
mission through a right angle junction of the airframe framing is 
an important case for consideration in performing such meas- 
urements . 

While the emphasis for the present study is on the coupling 
between structural subsystems, coupling loss factors between 
acoustic spaces and between structures and acoustic spaces are 
also readily evaluated experimentally. The experiments for the 
different types of subystems are, in general, the same. The 
source subsystem is excitated with a shaker or speaker and the 
space average response levels are measured in both the source and 
receiving systems. An independent measurement of the dissipation 
of the receiving subsystem is also required for the evaluation of 
the coupling loss factor. 

The basis for the experimental determination is described with 
reference to the energy balance equations for two subsystems where 
onlv one is excited. The power balance for the unexcited sub- 
system yields the following result. [17]; 
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N , N - « of modes in excitation 
^ ^ baodvidtli 


^ - coupling loss factor 

n - (lampiiiK loss factor in receiving 
^ subsystem 


Measured quantities include the total energies E ^oh, ^s. tot' 
the damping loss factor t| . The mode count is' deteriifineQ from 
appropriate analytical expressions for the particular subsystems. 
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The total energies are determined from the space average response 
level! Fira^lcoustic subsystem the total energy is the energy 
density, D, times the volume of the space, V: 


= VxD 

tot 


(132) 


where: 
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D = TT 
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<p2>^ space average mean square pressure 


(133) 


For a structural subsystem the total energy is related to the 
measured vibration response according to: 
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- surface density of structure at i^^ 
vibration measurement location 

- surface area associated with i^^ 
measurement location 
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measured sguare vibration level at 
i”" Location 


For a uniform structure with evenly spaced measurement locations 
this simplifies to: 


,,vi b 


M <V'2> 

S X 


(135) 


where : 


M - total mass of subsystem 
s 

<V2>^ - space average mean square velocity 
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The accuracy of using the above equations to experimentally 
determine the coupling loss factor is dependent on the damping 
level in the receiving subsystem relative to the coupling loss 
factor. If the damping, n is small compared to q /it is not 
feasible to solve Equation (131) for the couplin§' ^oss factor 
''r s' is ^ condition where the modal energies of the sub- 

systems are equal. It is often necessary in performing a labora- 
tory measurement according to Equation (131) to artificially 
increase the damping in the receiving subsystem to improve the 
measurement accuracy. For acoustic systems this involves the use 
of absorbing materials such as fiberglass or open cell foam. 
Damping treatments for structures typically involve added lossy 
viscoelastic materials in the form of free extensional or con- 
strained layers. 

An alternative approach for an experimental evaluation of the 
coupling loss factor between acoustic spaces is to measure the 
transmission loss of the intervening panel surfaces. The coupling 
loss factor is then determined from the measured transmission loss 
according to Eqs. (66) and (67). Three techniques are applicable 
for measuring the transmission loss: 1) the room acoustics 

method, 2) the near field pressure method, and 3) the acoustic 
intensity method. 

For all three techniques a speaker is placed in the source space 
to goner^te the noise field. The acoustic intensity incident on 
the intervening panel between the source and receiving space is 
evaluated for all three techniques according to the following: 

<p^^> 

^ i nc 4pc (136) 


“ sp.ice avt'raj'.e int'.ui S( 4 uare pressure in 
source space 


The techniques differ in the manner in which the transmitted 
acoustic intensity is evaluated. 

For the room acoustics method a prior calibration of the receiving 
space is performed which relates the measured space average mean 
square pressure levels to the power transmitted into the space 
under steady state conditions. The prior calibration involves a 
reverberation decay time measurement to quantify the energy 
absorption in the space. 
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The difference in measured space average pressure lev#»ic; in 


TL = NR + 10 Log,^ ^ 
*10 R 


(137) 


TL - transmission loss 
NR = noise reduction 

NR = <p|> - <p 2 > = difference in space average 

pressure levels in source (s) 
and receiving (r) spaces 

= area of intervening panel 
R = room constant for receiving space [19J 


A key factor in applying this method for helicopter spaces is 
It evaluates the total power transmitted into ?he JeceiSfng spacf 

intere^'s^ For° through the intervening oanel of 

i-h<^ the transmission loss values to accurately reflect 

m?ni»Sm transmission must £e 


estimations of ths transmitted 
^"tough the panel of interest and are not sensitive to 
flanking transmission errors. in the near field pressure method 
microphone is swept at close distance over the receiving side of 
near^f^eid transmitted intensity is related to the average 

preLion: pressure according to the following approximate el- 
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the^^oa^^^ and average near field pressure over the surface of 

radiation ohana/'^-Sl ^pineal factor that depends on the 

curvature of the surface including the effects of 

curvature. Typical values for are in the range from 2 to 4. 
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For aircraft structures additional empirical information may be 
required to specify the value of K^. 

The advent of small, high speed computers capable of efficiently 
performing digital signal processing operations has resulted in 
the implementation of methods for measuring acoustic intensity. A 
two microphone probe is used to measure the intensity of the 
radiating noise field in the vicinity of a vibrating surface. A 
more detailed discussion of the theoretical basis of the intensity 
probe is given in Appendix C. The application of the intensity 
method for measuring the transmission loss of aircraft panels is 
discussed in References (25, 26]. 

The intensity probe is used to measure the net radiated intensity 
from the intervening panel by scanning over the surface at a 
sufficient number of locations. The transmission loss is then 
determined by the ratio of transmitted to incident intensities. 
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TL = -10 Logio (T) 


(140) 


where I is evaluated from Eq. (136) and is the net power 
transmitSid through the panel as measurea^'wfth the intensity 
probe . 

The intensity techniques can also be used to evaluate coupling 
loss factors between panels and the adjacent acoustic spaces. 
These evaluations are based on an experimental determination of 
the radiation efficiency for the panel. With the panel excited by 
a noise field on the opposite side from the space of interest, the 
power radiated by the panel is related to the radiation effi- 
ciency, t>y the following: 

"rad = "rad 

where <V2> is the space average mean squared velocity level on the 
panel. ?he radiated power is determined by the acoustic intensity 
method accounting for the total radiating area of the panel. The 
experimentally evaluated radiation efficiency is then related to 
the coupling loss factor accordin<j to Eq. (68) as follows: 
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Airframe Vibration Transmission Measurements 

vibration transmission measurements were carried out 
on an S-76 airframe with the following objectives in mind: 

A) Familiarization with S-76 airframe to assist in design 
and implementation of full Phase II ground test measure- 
ment program. 

- Identify measurement requirements 

• source excitation 

• sensor location 

• data acquisition and reduction 

B) Provide limited quantitative descriptions of vibration 
transmission to assist in development of S-76 airframe 
model . 

Help identify important structural elements and 

paths . 

The measurements were performed on a partially assembled S-76 
airframe on the Sikorsky S-76 assembly line. All of the major 
structural components of the airframe were in place. The pilot/ 
copilot doors and windows had not yet been mounted. The gearbox 
was also not installed, and this conveniently allowed access for 
mounting a mechanical shaker to the gearbox attachment locations 
on the airframe. Control and electrical cables were not in- 
stalled. This allowed convenient access for mounting an acce- 
lerometer on the various frame and panel members of the airframe. 

For the measurements a mechanical shaker was mounted separately at 
two gearbox attachment locations. The shaker was excited with a 
broadband random noise signal in the frequency band from approxi- 
mately 200 - 7000 Hz. One lightweight accelerometer was per- 
manently mounted at the source location while a second accelero- 
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meter was moved to different locations on the frames and panels of 
the airframe. A total of 40 different locations were used for the 
measurements. The locations covered the structure around the 
cabin including overhead panels and frames, side panels, windows, 
doors, etc. and extended from near the source in the rear of the 
cabin forward to the nose structure. As a result of the symmetry 
of the airframe structure, the measurement locations were limited 
to the middle and left side of the cabin. The shaker positions 
included the rear most attachment point on the same side as the 
measurement locations and the forward point on the other side of 
the cabin. Lightweight .5 gm accelerometers were used so as not 
to influence the vibration levels being measured. 

The data were acquired and processed with a two-channel FFT 
analyzer which computed the transfer function between the vibra- 
tion levels at the source and measurement locations. Coherence 
function measurements also were recorded to indicate frequencies 
where potential contamination by extraneous noise sources occur- 
red. 

A sample plot of these measurements is shown in Figure 25 with the 
shaker located at the rear attachment point on the left side <RAL) 
of the gearbox and the measurement location at a random point on 
the overhead panel on the left side of the cabin near the rear 
bulkhead. The overhead panel is of a composite honeycomb cons- 
truction. Refer to Table 1 for frame and panel designations. 

The narrowband plot of the transfer function is characterized by 
peaks and dips that are indicative of the resonant structure of 
the panel response plus the fact that the result involves no 
spatial averaging over the area of the panel . The coherence 
levels are in general high, indicative of a large signal to noise 
ratio. The dips in the coherence function correspond to the dips 
in the transfer function. Dips occur where the accelerometer is 
located at a node of individual panel resonances. The measured 
panel signal at the dip is dominated by extraneous noise that is 
uncorrelated with the source signal. 

As shown in Figure 26 the overall coherence levels decrease for 
panels further from the source as a result of the decreasing 
vibration transmission. The transfer function magnitude is 
evaluated from the cross spectrum of the source and measured panel 
levels, and therefore, discriminates against noise in the panel 
levels. Even though the coherence, in general, is decreasing due 
to noise the accuracy of the transfer function estimate is main- 
tained. 
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The variation in overhead panel levels is shown in Figure 27 for 
the shaker at the front attachment point on the right side (FAR). 
The overhead panel levels were measured on the left side of the 
cabin. A comparison of levels on different panels is influenced 
by differences in the mechanical characteristics of the pani.1. 
The thicknesses of the composite honeycomb panels vary from 12 to 
25mm throughout the cabin. 

In an SEA model of vibration transmission, the modal energy is the 
important quantity in determining the flow of vibratory power. 
Differences in modal energies between attached structures result 
in a net flow of power. An equilibrium results when the modal 
energies are identical and therefore no net power flow exists. 
The modal energy is scaled by the space average vibration level 
where the constant of proportionality is dependent on panel 
parameters including thickness. The result of this is that a 
comparison of vibration levels is not sufficient for charac- 
terizing the vibration transmission when the parameters differ 
from panel to panel. 
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Figure 25. Typical Transfer and Coherence Function Measureiiient. 
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Figure 26. Variation of Coherence Function Levels for Overhead 
Panels vs. RAL Excitation. 
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Figure 27. Variation of Transfer Function Levels for Overhead 
Panels vs. FAR Excitation. 



From Figure 27 the variation in level from the rear of the cabin 
forward to locations at the copilot’s position is relatively 
minor, on the order of 5 - 7 dB. The absolute magnitudes of the 
transfer functions can be significantly greater than 1. This 
depends on the fact that the gearbox is attached to relatively 
heavy and stiff framing. If the framing is well coupled to the 
panels so that the modal energies are comparable, then the vibra- 
tion levels will be greater on the lighter weight panels. 

Differences in levels for individual panels for the excitations at 
two different gearbox attachment points (RAL vs. FAR) are shown in 
Figure 28. These differences are more pronounced for panels 
further aft closer to the source attachment points. For the 
forward panels the levels are more comparable. The differences 
are more accentuated at higher frequency. The source on the left 
side (where the panels are located) is closer and more directly 
coupled into the panels in comparison with the attachment location 
on the right side. 

Transfer functions for important panel surfaces in the rear of the 
cabin are shown in Figures 29 and 30. The drip pan is the over- 
head panel immediately beneath the gearbox that is supported from 
the framing to which the gearbox is directly mounted. The rear 
passenger bulkhead is also directly coupled to the same framing. 
Both panels are more equally excited by the two attachment point 
locations at lower frequency. At higher frequencies the rear 
attachment point is more effective than the forward location in 
exciting both the bulkhead and drip pan. 

The variations in vibration levels on the webs of the main fore/ 
aft frame members are shown in Figures 31 and 32. The rear 
attachment location (RAL), Figure 31, is at the aft-most point of 
the fore/aft frame. The levels are greatest on the first frame 
section, particularly at higher frequencies. The decrease in 
level forward of the gearbox is not significant. In fact there is 
an increase in level for the furthermost forward frame section. 
The decreasing weight of the frame members moving forward contri- 
butes to this result. 

The variation in vibration level for the right side source loca- 
tion (FAR), Figure 32, also does not show significant changes in 
level. The source is directly coupled into the cross frame 
between fore/aft frame sections one and two. Levels in frame 
section one are smaller than for section two because frame section 
one, which directly supports the gearbox, is of a heavier cons- 
truction. 
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Figure 28. 


Effect of Excitation Point on Panel Transfer Function. 
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Figure 35, 


Variation of Vertical Transfer Function Along the Ma 
Longitudinal Frame for FAR Excitation 
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Effect of Excitation Point on 
Transfer Function. 
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The flight test measurements are primarily intended to provide 
source level data for both acoustic and vibration sources. The 
measured source levels will be used to scale the inputs to the 

Recently developed partial coherence techniques 
will be applied to in-flight data to determine the contribution of 
each source to cabin levels. This will provide an independant 
check of the overall model prediction. 

An important additional objective of in-flight measurements is to 
detailed descriptions of the cabin noise environment for 
operating conditions. The descriptions involve 
detailed measurements of c^in pressure levels at the different 
passenger/crew locations including an evaluation of spacial 
level near the different passenger locations and in 
the vicinity of important radiating panels. It also includes 
extensive measurements on panels radiating directly into the cabin 
including overhead panels, the webs of main frames, door panels 
and windows, etc. Vibration and acoustic levels will be measured 
at intermediate locations on the airframe and in intermediate 
acoustic spaces, for example the luggage compartment. The mea- 
surements in the cabin and at intermediate locations provide 

of levels throughout the helicopter for comparison 
with the model predictions. 

Several flight operating conditions will be evaluated during the 
measurements, including: 

• 75 meters per second level flight 

• level flight at most efficient forward speed 

( longest range ) 

• hover 

Other details of the flight condition as regards weight, weight 

distribution, etc., will be consistent with standard flight test 

requirements. 

Initial survey measurements will be performed by moving a single 
microphone or accelerometer from location to location. A refer- 
ence location will be monitored to verify that the operating 
conditions have remained stable. This reference location will be 
monitored during the fully instrumented flight tests to verify the 
consistency of flight condition. Analysis of this initial survey 
data will be used to select specific in-flight measurement loca- 
tions for the detailed data acquisition required for the analytic 
model verification. 
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flight test measurements will involve a 
acoustic and vibration sensors located throughout 

recoiiir recordings on a 14 channel \ape 

recorder. Since the Phase III isolator design is likelv to 

involve materials with temperature sensitive characteristics 
temperature measurements near the isolator locations will he 

mated^for thl Wroximately 36 sensor locations is esti- 

mated for these flight test measurements. 
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Task 2. 0 flight test measurements . - These flight test measure- 
ments are primarily intended to determine source strength levels. 

2.1 Instriment model S-76 aircraft with accelerometers at key 
locations (gearbox, hydraulics, frames, windshield, etc.) 

2.2 Perform 1 hr. preliminary vibratory/acoustic survey to 
identify additional in-flight hot spots 

2.3 Acquire approximately 6 hours of flight data in various 
conditions measuring vibration and acoustic levels at ap- 
proximately 36 locations 

2.4 Acquire approximately 3 hours of flight data (subsequent to 
SEA model modification). Measurement locations will be 
chosen to verify model modifications and refinements. 


Model Comparisons and Refinements 

p.e primary objective of Phase II is the validation of the cabin 
interior noise model developed in Phase I . The development of a 
noise model for a system with the complexity of a helicopter 
necessitates an iterative process in order to establish confidence 
in the model. Measured data for a representative helicopter is 
compared with preliminary model predictions in order to identify 
weaknesses in the model which require additional refinements*. 

The validation process involves comparisons of the reduced flight 
and ground te«?t data with predictions from the model. The model 
is exercised for original estimates of parameter values in making 
an initial comparison with measured data. Subsequently, addi- 
tional predictions are generated to determine the sensitivity of 
the predicted results to variations in key parameter values. This 
also involves comparisons between predicted values for individual 
parameter values, example, coupling loss factors for specific 
structural connections with the results from specially designed 
laboratory experiments on structural subsections. 

The intended end result of the Phase 1 1 measurements and compari- 
sons is a cabin noise model that has been thoroughly validated and 
which IS readily applicable for evaluating the effectiveness of a 
variety of noise control treatments. 
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SUMMARY AND CONCLUDING REMARKS 
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sS££"£'is'""k 

^vss 
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under operational conditions, contribute importantly to one's past 
experience in initially developing the model and subsequently in 
refining it. The second model verification phase of the NASA 
program will provide important measurement inputs in demonstrating 
tne usefulness of a Statistical Energy Analysis approach for the 
complicated problem of modeling the cabin noise environment of 
helicopters . 
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Details of Source/Path Identification Techniques 

temp oral discrimination . - The most direct 
method of identifying transmission paths of a non-dispersive 
system is by exciting the system with impulses (using either an 
externally applied signal or internal impulsive sources witSn the 
arrival' monitors the response at a point and detects the 

the impulse. Since the wavespeed is constant, 
then a measurement of the travel time leads to an estimate of the 
propagation distance between the points if the wavespeed is known 

betw^er^^^ir^noinr**^^® wavespeed if the spatial pith separatioA 
between the points is known. Examination of later impulses 

information on the first few reflections within the 
paths can be identified and their relative 

amplitudes measured. 

Problems arise with all time delay methods when the arrival times 
for successive propagation paths become closely spaced relative to 

For vibration transmissiojl L I 
helicopter airframe this may constitute a serious limitation 
because the structure is relatively compact and path l2Sg?j; 

small. The technique would involve excitating the 
gearbox attachment points with tone bursts or impulses and moni- 

acoustic signals at different passenger 
locations. Vibration levels on different panel surfaces of the 
cabin interior may also be monitored. 

— correlatio n analysis - A second method 
of temporal discrimination involves the use of correlation analy- 

anw function is computed between the 

t-h t response point if the source location is known, or 
between the two response points. If the source excitation is 

band-limited noise, then the cross-correlation 
^ maximum at a delay time which is equal to the time 

imio^ta^nt 5®^'^®en the two response points. The most 

thA between the time delay measurement method and 

ncff.n analysis method is that the time delay technique 

is pulse excitation, whereas the correlation technique 

ov steady-state, random excitation. if the source 

excitation is broad-band and produced within the system, then 
correlation analysis is often the more appropriate"^ technique. 

Whose path reftonses x(t) and y,t.x, 


0 


(143) 


Rxy(i) = XY cos[n(T-d/c)(fj^+fy)] 


sin(n(T-d/c)(fy-fj^) ] 
n(T-d/c)(fy-fj^) 


where X and Y are the signal amplitudes. The cross-correlation 
function is a maximum at x=d/c, the propagation time between the 
two response points, and the envelope of the peak is more sharply 
defined as the noise bandwidth (f^-f.) increases. Therefore, 
determination of the peaks in the Spatial cross-correlation 
function enable a measurement of the time delay between the 
signals . 

Now suppose we have a system with N paths such that the path 
length d- (i = 1, 2, N) is different for each path. The 

resulting cross-correlation function between the input and output 
has N peaks as shown in Figure 37. The peaks in the cross-cor- 
relation function have exactly the same form as in Equation (143), 
so that the delay time for each path corresponds to the time at 
which each peak of the cross-correlation function is located. The 
amplitudes of each peak depend upon the relative magnitudes of the 
transmission coefficients of the various paths. 

Figure 37 illustrates that various transmission paths can be 
identified only if the delay times between any pair of paths are 
sufficiently different that the peaks can be distinguished. Let 
Ti and i 2 be the delay times for two of the transmission paths 
between the two response points. The criterion [27] such that 
these paths may be distinguished using correlation analysis is 
that the difference in these delay times satisfy the following 
inequality: 

1^2 - ti I (fy - fj^) > 1 (144) 

The normalized correlation coefficient C (x) is defined to be, 

xy 

c (X) = R (t)/[R (0) R (o)l^ (145) 

xy xy XX yy 

where R is defined in Equation (143) and R and R are the 
input attd output autocorrelation functions.* The correlation 
coefficient for the time delay between transmission paths is given 
by Equation ( 143 ) with the signal magnitudes X and Y set equal to 
unity. Thus, the correlation coefficient evaluated at the delay 
time is equal to 

C,„(x=d/c) = 1 (146) 

xy 
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Cross corrdalion in multiple-path sysic 


Figure 37. Cross Correlation in Multiple-Path Systeii 



It is important to note that the amplitude of the correlation 
coefficient is independent of the transmission path characteris- 
tics for non-dispersive systems. This is not the case for disper- 
sive systems. Figure 38 shows the correlation coefficient for a 
dispersive system. The time width At between the adjacent zeros 
of the correlation coefficient on either side of the central peak 
is illustrated on the figure, and this time width is equal to 


At = i/(fy - = i/Af (147) 

where Af is the frequency bandwidth of the excitation. The cor- 
relation coefficient peak is most clearly resolved by increasing 
the frequency bandwidth of the excitation. Equation (147) also 
helps illustrate why the condition expressed in Equation (144) 
must be satisfied for multiple transmission path systems in order 
to resolve individual transmission paths. 

Temporal discrimination: the frequency domain . - The time delay 
between various transmission paths also may be measured in the 
frequency domain. A common method is to use the cross-spectral 
density function. Consider a signal s impinging on each of two 
response points with excitations x(t) and y(t) such that a noise 
signal is present at each response point as follows: 


x(t) = s(t) + m(t) (148) 

y(t) = s(t+t) + n(t) (149) 


There is a time delay t of the signal at y relative to x. The 

noise signals m(t) and n(t) are presumed to be uncorrelated with 

s(t) and with respect to each other. The cross-spectral density 

G (f) between x and y is 
xy 

0,/f) = C^^(f) (150) 
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Fi(ju»e 3o. Cross Correlation of Dis()ersive Waves. 



The autospectral (power) density of signal s is G 
real non-negative quantity. Therefore, the phase®!' 

xy 


and it is a 

of Sy 


4> 

xy 


U»T 


(151) 


between the transmission paths is simply obtained 
^ency *** araph of phase versus L|ular i«f 

It is soiuetimes possible to improve the -4 i 


Y (f) = 
*xy 


G 


XX yy 


(152) 


coLreTce" djs^issed'”'’*^^* 

function for the system defined by Equations (148? andTl49?^is 



G e 

J{G + G 
ss inm 


jUJT 

m 


ss 


Tin 


(153) 


of ”°^®® autospectral densities. The phase 

iibherence function is precisely equal to the phase 
or^JUf^ Equation (151), so that the complex coherence function 

oSlurs*\!hen^toi"‘;\^rre ‘h® f«o frequency domain methods 

Sign^ s Tte sTonai ore partially correlated with the 

Ta. Signal s in practice may not be a white noise 

well as a few^'strona^*<f^^^^'- ® ®^ 9 nal with noise components as 
nof Kaf^ea-^a- ■ sinusoidal Components. Such a signal will 

thl nLi independent relative to a noise signll Thus 

the phase of such a cross-spectral density is not proportional to 
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relatlo^n^o “th 4 the'" "*013^ soM™^• noise*'*'i 

Shilh"is norSI“!n?«Lt^"''Thl"cr^ 

the complex coherence 

response. The complex coS«eLe wtiSS 'is thf ratio*®^^^^ 

aJ?ll;ell?a“"'ae^^??ie^^ T, Si °"^%roaSct °Jf Si 

:ScT-h-'trfi 5 

ssi. “r—rHEP H‘js 


xy 


(T) =-00 / “ W(f^ y (f) e‘2^j^^df 


(154) 


hSSIeS=\-„r=hSoVa SS^noiSe SSnaSl%^roSuc^err 

orth'e'sYnuso'idT S" lddSSI"SV>| nolst"\e''adS"t ' 

iSsi 

L“ — -~t SI 

complex signals °f delay time for certain types of 
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•€lAHV€ IiMf OCIAY 
Croii correlatica of x(r) and y(0. 



tflATiVC TiMf OClAY 
N'orm.ilixed SCOT of x(/) and y(0. 


Figure 39 Normalized SCOT. 



ot transmission paths in dispersive lystems such^i wails 

in structures is complicated by the Lcf that the 

the wave excitation is not constant with frequency ^ As a ?llult^ 

there re not a well-defined delay tire between ?^ISs™on patos 

prevToreirSiscusred tl°" 

‘SS ‘ ~:c™ i‘:s 

Cross-correlati on analysis for a d isDersivo svstem - Th«»r-^a =.v 
significant differences between the cross-co/relation 
dispersive versus non-dispersive sy«IL The °S 

r1:«L^rrhrt\re"r?'^‘‘°^ rr 'r«rly "rtoS'^ri .ll 

a sufficient effects can be minimized by choosing 
velocity is' nearly ^ions^anT^’^uc^t'^^ 

time width \ systems notes that the 

inversely proportionaf\o °the ^requenc^^b^^^^^^ 

tL timl ^r which ^ estimated by measurement of 

exi??atrr^\F°”- 

ff ?f Jn^rrotilir^L "f^riurn^ ^Cn^-rth "°to^Telrer Si 

delay time. increases the uncertainty in estimating the 

This problem becomes more pronounced in a Qv«t*»m e.,^v, u i 

iim'issSrvk£7;^s 
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1 


lower and upper cutoff^fre<JiencieW ^and 

a-s^n\Te^TrSi3sL^^ 

(dispersive) bending waves on an lircra^ft fuc Propagation of 
excitation point to a distant noint on^^K from a source 

correlation coefficient c cross- 

cross-correlation function^^is ' ^ helo}>i as the normalized 

R ( 0 )]^ 


(1/Af) ^ / cos[27if(T - d/c (f))] df 
jL p 


(155) 


frequency bandwidthT^and*c°’( £)’*'is difference, af is the 

which varies with frequency? ** dispersive phase velocity 

T^^pCse'vlUci?? f«"a‘be®nd\i!f’'®‘* bending waves. 

and longitudinal 7 hasf velocity c" ir" thilkness h 

L 


Cp(f) = [uihCj^/^12]^ 


(156) 


Define the center frequencv f of fho 

^equency of the frequency band as follows: 


^0 = ^ 


(157) 


Cl • V - / 

choice of phase velnn*i t\r *i-k^ 

cient in Equation (15^5) is found to' b? '"^“ss-cofrelation coeffi- 


<.,/t) = |2y(Afi,^)| (cos|fi(Ci(z^) - CI(Zj^)| 

+ sin(fi(Si(Zy) - S:(Zj^))} 


(158) 


0 


where 


^ ^ - 3/4)2 + 3/8J 

^ ^ ■ 3/4] 


(159) 

(160) 


^h© t©i*in z is evaluat'f^H ai- 

exes in Equation (158). The Fresn^^^®^ frequen- 

[30] are Ci(z) and Si (z), respecUv^^ integrals 

evaluation ?f^?jrFre^snel^inTeV^^^^ '^The'^D^ ^ detailed 

relation function (or cross-coSl ^ i cross-cor- 

txme Xg equal to the group delay dehned a'^^fouO^ws''^^ ^ 


The group delOy^^is*^^thl Ume^^iOte^val^ energy propagates, 
distance d. The group velocitv iO^O^^ 1 ^ energy to travel e 
in a dispersive system; it is given by phase velocity 


C»UI/<7K 


F H- 

giJeniriqiluli®*"?} anr?he"g?oJr^“o?itris'’*’^^® velocity is 


C 

8 
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non-dispersive syLeroeSjs at f Peak for a 

phase delay time is equal to the arouo^doi^®^ 

group velocities are equal in a Ln®-Tlpersfve syTtem’;'^ 

Iht g'«ss"cot«u\iol <« wall as 

given by «ir-cient) for a dispersive system is 




_ HU- (164) 

The center frequency of the exc-i •^;4•t-•i r^n •; r 

tp IS defined in terms ofthf phase vainai'^v ^ P*'®*® 

difference d by, pnase velocity c and path lengt 


(165) 


The width of the time window in Equation (164) has been defined to 
be the time interval in which the cross-correlation coefficient in 
Equation (158) is greater than or equal to one-half of its peak 
value. There are two differences between the width of the time 
window for dispersive and non-dispersive systems. The width of 
the cross-correlation peak for a non-dispersive system is indepen- 
dent of the delay time whereas for the dispersive system it 
depends on the phase delay time. Secondly, the peak width of the 
cross-correlation function is more spread out in a dispersive 
system than a non-dispersive system. 

The peak amplitude of the cross-correlation coefficient for a 
non-dispersive system is equal to unity independent of time delay 
fr6^u®ncy bandwidth. This is not the case for dispersive 
systems. The mean amplitude of the cross-correlation coefficient 
for a dispersive system within the window width corresponding to a 
single peak is 


c 

xy 


= 2 (Vp) 


- 1/2 


(f()/Af) 


( 166 ) 


This demonstrates that the cross-correlation coefficient can have 
a value much larger than unity, particularly for narrow bandwidth 
signals. Equation (166) is valid when the path length difference 
is large relative to a flexural wavelength. As the path length 
difference becomes small relative to a wavelength the cross- 
correlation peak value approaches unity. 

The results of this section are summarized in the statement that 
the cross-correlation coefficient contains information on the 
dispersive system transmission path characteristics with respect 
to the time at which the peak occurs, the time width of the peak, 
and the mean amplitude of the peak. The cross-correlation coeffi- 
cient contains information on a non-dispersive system only with 
respect to the time at which an individual peak occurs. 

Cepstral analysis . - This section describes a nonlinear filtering 
technique known as cepstral analysis. Due to the advent of small 
computers with extensive computational power, practical applica- 
tions of cepstral analysis for filtering and signal recovery in 
the laboratory have become possible. Cepstral techniques have 
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over tlie more traditional linear signal processing 
techniques in the areas of dereverberation or spectral smoothing, 
the analysis of periodic signals, source identification, and 
transmission path analysis. 

Cepstral analysis may be used to detect and/or remove echoes from 
a signal, and can thus be employed to compute the free-field 
response from the measured reverberant field response. Of parti- 
cular interest as regards helicopter applications are the capabi- 
lities of cepstral analysis for periodic signals. The periodic 
components of the spectrum produced by harmonics of a particular 
gearbox rotation rate is one such periodic signal. All of the 
harmonics corresponding to a particular rotation rate are mapped 
by the cepstrum into a single peak. Filtering in the cepstral 
domain thus isolates the source characteristics of the various 
gearbox sources. The cepstrum may under certain conditions be 
used to separate source and transmission path effects. This is 
because the cepstrum converts the convolution of source/path 
characteristics in the time domain into simple addition in the 
cepstral domain. 

Future cepstral applications to systems, such as a helicopter 
frame exhibiting both dispersion and reverberation, depend upon 
the development of more powerful techniques than exist at present. 
If the dispersive characteristics of the structure can be ade- 
quately modeled, then a combination of cepstral techniques with 
inverse filtering (used to remove dispersive effects) can be 
employed to study structural transmission paths. 

The complex cepstrum k (q) of a .signal y(t), whose Fourier trans- 
form IS Y(f), is the iriverse Fourier transform of the logarithm of 
Y(f). Since the transform Y(f) is a complex quantity, then the 
appropriate logarithm is the complex logarithm. 

The transform Y(f) is written in terms of its magnitude and phase 
as follows: 


Y(f) = |Y(f)| 

and the complex logarithm of Y(f) is 
log Y(f ) = log |Y(f ) I + j({) 


(167) 


(168) 


The real part of the above involves the logarithm of a real 
nur^er, and this is well-defined. The imaginary part is not 
uniquely defined since addition of 2n radians to the phase of Y( f ) 


I 
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necessary 

to restrict the range of the phase to 2/i. Under this condition 
the complex cepstrxm Jc (q) is defined to be the inverse Fourier 
transform of the logari^^hm of Y(f) as follows: Fourier 


= S log Y(f) df 


( 169 ) 


The parameter q plays a role analogous to that of ^-ime in the 

^eSstrui""?? ti^^^sform. During the corrnutation of Se 

cepstrum It is necessary to unwrap the phase; i.e., the phase must 

nutation function of frequency. Information cn tare com- 

putation of the complex cepstrum is foun^ in the references. 

The power cepstrum c (q) is defined to be the inverse Fourier 

transform of the log^ithm of the autospectral density G (f) of 
the signal Y(f) as follows: aensity t-yyU) of 


lo8(Gyy(f)l 


( 170 ) 


autospectral density is a real, non-negative function, 
then the logarithm in Equation (170) is the real logarithm. 

are real quantities. This is obvious in the case of the power 

complex cepstrum is real depenIs 
upon the symmetry properties of Equation (169) [31] • 

The complex cepstrum is more general since the original time 
domain signal may be completely reconstructed, whereas this 
reconstruction is not possible with the power cepstrum. 

Cepstral analysis has usefully been applied to signal deconvo- 

The^ input and^ input, single output Unear system, 

ine input and output time domain signals x(t) and y(t), respec- 
tively, are related by the following convolution: ^ ^ 
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y(t) - h(t) ■' x(t) 


(171) 


y(t) =_J' h(T) x(t-T) dT 


(172) 


(171) is a short-hand notation for Equation (172). The 
impulse response function of the system is h(t). 

domain Lecomes addition in the cepstral 
omain. It is this property which makes cepstral analysis useful. 

the source the cepstral domain to focus attention on either 

(x(t)) or transmission path (h(t)) characteristics. A 
periodic source excitation produces a cepstrum with a peak at a 

nVfll value of q which is directly related to the repetition rate 
of the time domain signal. 

filtering techrdques in the cepstral 
rhaf frequency domain. For example, assume 

input x(t) to the single input, single output system is a 

also produced by gearbox noise. Assume 

also that the impulse response function of the system is slowly 

^®P®tition rate of the pulse. The convo- 
^ impulse response function produces a 

5 ”^*^se high q values are predominantly due to the pulse 

values are due to the transfer function. Appli^ 

trai^ff^r cepstral domain isolates the 

transfer function contribution. 

domain signal x*(t), after filtering in the cepstral 
domain, may be recovered by use of the inverse cepstrum. Com- 

details^o/th^^ffi?^^- Phase considerations and the 

details of the filtering in the cepstral domain. 

isolate transmission paths by 
flni si^al from the first few reflections, or 

from the reverberant field as a whole. The main property of the 

^ '^hich allows this is that a periodic time signal produces 

a cepstrum with a single peak. This method has been applied f321 

radiated sound from a jet engine, since 
that study was interested in estimates of radiated power rather 

was\s^ed°”^^^^^^^°” temporal waveform, the power cepstrum 
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Figure 40a illustrates the radiated power spectrum, and the 
ripples in this spectrum are due to ground reflection. The power 
cepstrum corresponding to this signal is in Figure 40b, and the 
ground reflection is represented by the very strong peak in the 
cepstrum. This peak has been filtered out of the cepstrum in 
Figure 40c, and the resulting inverse cepstrum leads to the 
modified power spectrum shown in Figure 40d. This power spectrum 
is compared to the power spectrum measured under anechoic condi- 
tions (i.e., when measured in a special test facility which 
virtually eliminates all reflections in the bandwidth of in- 
terest), and the agreement is excellent, Syed et al. [32] remark 
that the cepstral technique is applicable even in the presence of 
wind, turbulence, or temperature gradients near the ground. They 
found it useful in both narrow-band and third octave averaging. 
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APPENDIX B 


Frame Junction Transmission Model 


part of the^overall^odel^of^tL^i^^^ k” integral 

The fundamental feature^ o/ Se a^rfr^e 

transmission through the fr^mir^rr is the vibration 

structure and the skeletf[ including the gearbox support 

Sikorsky s-76, the f ?^inf cons?sTs''°o ? tJe 

with wrap-around tran^Jerle ?rames ^ frames 

and side panels of the cabin Thi ^ extend between the door 
tical model describing the vibrI?tL ^he analy- 

sections joined at right angles. ” transmission between frame 

° s ?h°or? n ngurff"\f 

are treated as beams that deform in h^nH • ^^^es 

perpendicular to the frame axis and ?n ^ ^ two axes 

axis. For both bending and tor<5inn.i about the frame 

section does not deform. At Mah^r f ™°^ions the frame cross- 
deformation occurs and the dvnAmiV? k k cross-sectional 

of flat plates. Longit^diniTTotion ‘hat 

s??ffe";\r„%-d‘- ort^eV 

contribute significantly to the overa^f VilrTuorirSJ^*?o^° 

is treated as a'*"poi'nt comiection^”'!t ^hown in Figure 41. it 
transmission through the junction frequency. Vibration 

motions. The exclusion 0 ^ 1000 ! a result of junction 
restricts the allowed junction motion from txhe model 

about the three coordinate axes Ind rotations 

perpendicular to the plane formed bv f ^J^anslational motion 

translations are suppressed bv th#» i^y rames. The in-plane 
longitudinal motion in Jh^ frames m-plane impedance for 

which can^occu*’^ between fr^es°*' motions illustrates the coupling 
the 2 axis inducefbeSiS^rformatfon 

aligned perpendicular to the 2 axis and i-*' ^fames 1 and 3 that are 

2 and 4 which lie along tbt I axis 1, ‘■^^mes 

in the Y direction induces LndlJo translation 

tached frames. ■•nauces Pending deformation in all four at- 
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Junction Displacements 



out of plane translation 
rotations about the coordinate axes 


Figure 41 


Coo^dLaL^S^tem“*'^^^°'" Showing 





Incident bending motion in frame 1 involves both vertical motions 
and rotational motions about the Z axis. The induced junction 
motions couple the incident bending motion directly into bending 
motion in frames 1 and 3 and torsion and bending motions in frames 
2 and 4. In addition, if frames 2 and 4 are not identical, 
junction asymmetries result in a coupling of the out of plane 
junction motion into a rotation about the X axis and torsional 
motions of frames 1 and 3. The coupling due to junction asymmetry 
is the result of unbalanced reaction moments in frames on opposite 
sides of the junction. 

Incident torsional motion in frame 1 couples directly into torsion 
of frame 3 and bending of frames 2 and 4. Junction asymmetries 
com,^lete the coupling of the incident torsional motion into all 
the possible transmitted motions, as is generally the case. 

In-plane bending about the Y axis couples only into in-plane 
bending of the attached frames. 

Transmission coefficient description of junction dynamics . — The 
SEA model quantifies the coupling considerations described above 
in terms of a set of coupling loss factors. A convenient way of 
analytically modeling the coupling loss factors is in terms of 
energy transmission coefficients between the incident and trans- 
mitted motions. The transmission coefficient is defined as the 
ratio of transmitted to incident powers for the types of motion 
under consideration. 

The transmission coefficient is analytically evaluated for the 
case where the attached frames are infinitely extended in length. 
The dynamic behavior at the junction is the same, but this simple 
formalism removes the need to consider reflected energy from other 
junctions in the attached frames. Multiply reflected energy 
leads to the build up of a reverberant vibrational energy field 
which is accounted for in the SEA model in the difference between 
the coupling loss factor and transmission coefficient. 

As shown in Figure 41, incident bending or torsional energy propa- 
gates towards the junction where it is either transmitted or 
reflected. The solutions for the transmitted and reflected 
energies are obtained by considering the superposition of blocked 
and radiated cases as shown in Table 5. In the blocked case the 
junction is prohibited from undergoing both rotational and trans- 
lational motions. The blocked forces and moments, which would be 
required to balance the effects of the incident wave and prohibit 
any motion at the junction, are given in Table 6 for both bending 
and torsional wave incidence. By definition in the blocked case, 
no vibration is transmitted to the other frames. 
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junction motions equal thSle ^for removed and the 
junction motions result in k full problem. The non-zero 

attached frames. For bendiIJrr ^ torsional waves in the 

is localized to the rS “»otion oc?ur? onl 

tially while the other propagates unatte^u a?' decays exponen- 
K the junction absence of 

which IS transmitted into the f^’ame^^ latter determines the power 

rype of motion*ejiate^trie junction for each frame and 

at the junction. For bendin^mo^Lr^J.!”^ to forces and moments 
forces and moments to transveree^e^ri^^ impedances relate shear 
For torsional motions the relalLr^^ angular rotations, 

angular rotations. Analytic So?eSio only moments and 

given in Table 7. ^ expressions for the impedances are 

summing^'^the m'lJlJ^nS' and ®f^?cls “d^e^fiiTed^ b appropriately 

for each attached frame: efined by the impedance matrices 
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junction impedance matrix: 
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zdj 

m,n 

- 

z( 3 } 

m,n 
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Z^ 

m n 

X, 'y 

z(2i 

m,n 

- 

z(Ai 

m,n 
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where the superscript denotes the frame number. 


si^etric where the main diagonal is the sum of 

snu;r.« “ s“5‘ ■ 

g;£..-=-ss. s™ •;! 


incident motion " junction motions in terms 'of tlie 

infu^Lmfssrrcoiyfici'^r 'to ''refat"f^^h2^ *"-l“-ting 
ralirt°:rin“l!'‘='' i^utoVe*^\\\ct"d K." ”T°h1°p"|„e‘? 

?o"?io“ 


Bending motion 

Power associated with shear forces/transverse velocity 
Pf = 5 Re (F n*) 

Power associated with moment/angular velocity 
= I Re (M e*) 


(181) 


(182) 


real%Irt^”°^^® complex conjugate and Re refers to the 

Sfs.r’sr"a.T„';‘s:s s “ .r,S5.‘T“;s 

a power is the sum of terms due to forces and moments: 


abend 
■ total 


2 Re (2^) In^l + i Rc(2jj^) |9|^ + 1 


2 ^ ne ) 


(183) 

whLe°the'"''torsfo^,i^^^^^ expression involves only the moment term 

angula“;elo?r?r?s',iXl.riabu'l""''® 


.tors 


total 2 ^^t^ |8|^ 


(184) 


tAg=f«:e7s%Yven «”tlon'‘o* 


e = ik C 

inc inc 


(185) 


This relationship between the 
amplitudes yields the following 
on the junction: 


angular and transverse velocity' 
expression for the power incident 
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(186) 


'’inc * 5 + Ae (zjl (186) 

Sample evalu ation of transmission coefficient . - To illustrate the 
process of evaluating the transmission coefficient consider the 
case of an incident bending wave in frame 1. The bending involves 
rotations about the Z axis and transverse motions out of the plane 
?o simplify the algebra opposing frames (1,3) and 
( ,4) are identical so that off diagonal terms in the junction 
impedance matrix are zero. 

Summing the moments/forces due to the junction motions with the 
blocked forces/moments and setting the sums egual to zero yields 
the following: 


m i 

Zj inc 
m 

z 


(187) 




(188) 


0=0 

X 


(189) 


M : 0=0 

y y 


(190) 


The incident power is given by: 


‘’inc = 2 * Ae(2<'>)} 


(191) 


The power transmitted to frame 2 is obtained from the following 
equations for bending and torsional motion, .-espectively : 


I Re(z[2)) Ip |2 = 


|Re(zJ^^ icjh ^ 

i 2 

izj I 


(192) 


.2 

I 


(193) 


>( 2 ) 

tors 




|c“ 


m 


j |2 


izl I 


\L 


inc 


m 


The transmission coefficients between bending in frame 1 and 
bending and torsion in frame 2 are determined from the ratios of 
transmitted to incident powers; 


^(1,2) 

bend, bend 


Re(Z^^^)|cJ^|^ 
(Re(Z^^^) + k^ReCZ^^)} \Z^ 

z 


and 


^(1,2) 

bend , tors 


Re(z[^^|C^^|2 


{Re(Z^^^) + k^ReCZ^^'')} IZ^ 
t mm 


(194) 


(195) 


With the frame impedance expressions from Table 7 and the blocked 
force/moment expressions from Table 6, these become: 


4 Re(Z^^h Re(Z^^h 


IzJ 1^ 
y 


and 

,(1,2) 1. ReCzf >) Re(2<‘>) 

bend, tors 2 

m 

z 


^(1,2) 

bend , bend 


(196) 


(197) 


These expressions have a relatively simple form in terms of 
impedances for the individual frames. It should be noted that the 
expressions involve only moment or force impedance and not com~ 
binations of both. 

The final step involves evaluating the coupling loss factor in 
terms of the transmission coefficient. The coupling loss factor 
is 
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( 198 ) 


/ 1 , 2 ) 

bend , bend 


ReCZ^^b Re(z}'^^) 


/ 1 , 2 ) 
bend, tors 


Re(Zp^) Re(Z^^^) 

t Ijn ' 


IZ^ 


(199) 


where - mode density of frame 1 in bending. 

coupling 

coupling loss^?acto?s be^^eei obtained for thi 

IS symmetric and the frames identical. because the junction 

f^ame^r^i m^^re^^oVlTc^ii^^^ ^fh^ Sa^Lr^tt "d" 

bending motion for the symmetric j unctTon casr®"" 
junction motions, n and 6 contribute to Tkok^ ® non-zero 
propagates away frdm the function ^nto^ f bending energy which 
power now inclu'des terS dL to both" ° anT"e" transmitted 

moLnt‘’anrfMcl°LpeLnler"'^?hi°no“®*^‘^ contain both 

"cont^r^^ g 

tra’n^rS^ref/^^^^^^^^^ buj“th\“f“ " doterminaWo^Jf 

simply presented in te™s"'‘o\""iramf im^eL^^s ! 

£rames''®al*so‘' resS^t^ "on-identical 

transmission coefficients in expressions for the 

couples into motions th^t it d^dno^ incident motion now 

casL These lsSaIeS‘’bv Jhi° =y™>btric 

though the algebra is more involved expressions given, 
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Table 6. Blocked Force/Moment Expressions 


Incident Bending Wave: ^inc® 



Incident Torsional Wave: 0 + wt) 

X, me 



See Tables 7 and 8 for definition of parameters. 



ORIGINAL PAG;;: ::j 
OF POOR QUAirnr 

Table 7. Impedance Expressions for 


Bending Motion 



'v -I, 


Torsional and Bending Motions 

^ 2 ' ®2 ” rooTOent, angular 
rotation about z 
axis 


^y' Hy “ shear force, 
traverse 
displacement 


‘'f 

^m,r| 


m, n 


ro 


dot above symbol denotes 
time derivative 


7 - Bk 3 ^ 

- V 0 - 


"m,f^ 


Bk^ 

ut 


V _ Bk . 

■ w + i) 


Torsional Motion 

y 


— X 


M = 2 e 


s 


M 


L. 


inoynent about axis 
angular rotation about axis 


Zt = Pl^C 


P t 


see Table 8 for definrtion of parameters. 

S00 f 2 l If ^ 

evaluating impedances of mechlnical° s% terns procedures for 
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Table 8. Dynamics Parameters for Torsion and Bending of Frame 


Bending Wave Motion 

FT ^ 


k = w/C^ 


B = El 


“ bending wave speed 

- bending rigidity 

- bending wave number 
“ radian frequency 

- Young’s modulus 

- density 

- cross-sectional area 

- moment of inertia of 
cross-section about 
bending axis 


Torsional Wave Motion 


c = (~) 

t ’'Dl 


^t ~ torsional wave speed 
G - shear modulus 

J • torsional stiffness of 
cross-section 

Ip - polar moment of inertia 
of cross-section 

p • density 



0 
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friHbmfg^rc.i'^ghax^cJr^st^crwer^ 

function representative of the overh^aH f right angle frame 

S-76. The framina is a overhead framing on the Sikorsky 

T section flanges ^coMected bv wi»h°^?^^'^^^^°” ^pper and lower 

approximately I m Tnd the widto 0^^^^, ’'‘'® "^iSht is 

.05 m. Plate thicknesses varv fr-r?r5 flanges is approximately 

The calculated fra^e inLtS tZ- W^^^^^^tely .8 to 1.5 mm 
given in Table 9. For sTmpAc^^^^^ 

were taken to be identical^ frames at the junction 

tion reduces the number of non 7 t^rr> sy^etry at the junc- 

coupling loss factors. n-zero transmission coefficients and 

c“y'betw«y“nc^^^^^^ in Figure 42 for the 

source frame and bending and torsioMl*moy"'^^"® ">otron in the 

fiLnATy^g^Lte^^thSn"^"^^^^^ is signi- 

sional motions of the riaht^anoA^f either bending or tor- 

to the large bending rigiditv of^th^^f"*®® ’ is related 

of the overhead f rfmi if is of or bending rigidity 

in supporting the gearbox "and cabin ^ importance 

for out-of-plane bendina iIotioT frames are quite stiff 

motion. The coupling a?rosf r torsional 

because the cross frame impedance lunction is large 

mi??if "" -Pedance discontinuity for thfASf de^f "ieidtng 

to the 


' 1,2 


c 

2wL 


‘ 1.2 


( 200 ) 
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C = 


g vciwuxuy lor incidental motion 


L = 


2 Cjj for bending (Cj^-bending wavespeed) 
for torsion 

length of source frame 
w = frequency (rad/sec) 

cient values^i/^n^fA transmission coeffi- 

frequency dependencies in Figure 43 are thiiestir The different 
and the frequency dependence" of"c/fn EqLuoiuoo,"*'* 

nlw Cui“f ChCiourcf °fram"e j" .^““hi tying the po.er 

damping loss factor. The signifila^CoCth"" es che 

values can be understooH ^ ^ coupling loss factor 

factors. For riveted struc^urf^^ typical damping loss 

.01 are common a thT^v?! 

small in comparison with thf is 

motion into the in-line f^Le factor for bending 

because the enerqy which trancsmiVe ?°''’P^^iaon is not exact 

to reflection a^ revCrbCritlCn in CL'""” <*>'® 

accounted for in this calculation frame is not 

in solving the SEA matrix equa\?ons \Cr t\“"fuh^?st%TCnCr^ie1: 

anrtor\^"n^\""mo\°^?ns^Cre°"sU^^^ 

"motr ^isl^Siriy™d Cn"“"l, ^ ® >-h^fn. YnlZul 


187 


Parameter Values 


Table 9 Representative s-76 Frame 
Material - Aluminum 

E = 7.2 X IQio n/m2 
G = 2.7 X 10»o n/m2 
p = 2700 kg/m2 
Frame Geometry 



Parameter Values 

Moments of Inertia of the cross-section 
Ijj = 3.2 X lO'Sm^ 
ly = 4.29 X lO-^m-* 

Torsional stiffness of the cross-section 
J = 1.69 X lO’^m** 

Cross-Sectional Area 

= 4.43 X lO-^nji 

ngth of source frame between cross frame locations 2 and 3 



'iftimMilliiiiilii la tfilY'i III Ikimi ii 


POOR QU,,U,y 


Incident Motion: Bending, out-of-plane 


IN LINE FRAME: BENDING 


RIGHT ANGLE FRAME : BENDING 




RIGHT ANGLE FRAME: TORSION 



4 5 6 


Frcquonc') 


Figure 42. S-76 Frame Junction: Transmission Coefficients for 

Out-of-Plane Bending. 
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APPENDIX C 


Acoustic Intensity Method 

The acoustic intensity method involves scanning over a vibrating 
and radiating surface with a transducer which generates outputs 
that can be processed to give the acoustic intensity or power. 
The net acoustic power from the surface is obtained from an 
average of the intensity evaluated at close distances to the 
surface and at a sufficient number of locations over the radiating 
area. This method does not generate information characterizing 
the directivity of the acoustic power. 

It is important to distinguish between the intensity and pressure 
or velocity variables which are used as inputs and outputs. The 
intensity is a flow variable describing a net flow of energy 
through a surface in space. It is obtained as a time average of 
the product of in-phase components of pressure and velocity on the 
surface. For different surfaces which have coherent components of 
velocity an additional factor is involved in assessing the in- 
fluence of that coherence on the radiated intensity. This factor 
is that there is a mutual coupling between the radiation imped- 
ances of the two surfaces. That is, the pressure at surface one 
is influenced by the motion of surface one, but it is also in- 
fluenced by the motion of surface two as well as any other vibrat- 
ing surfaces. 

The definition of acoustic intensity in terms of acoustic pressure 
and flow is now introduced. The acoustic pressure p and particle 
velocity u in the r direction are related by 


= - O/P) /dt 8p/8r 


( 201 ) 


where the equilibrium density of air is p and the integration is 
over time. Thus measurement of the pressure gradient 3p/ar in the 
r direction enables a computation of the flow in that direction. 
The acoustic intensity I in the r direction is defined to be the 
time-average of the prodifct of pressure and flow as follows: 


I = <pu > 
r r 


( 202 ) 


where the brackets '<...>' denote a time average . 


Acoustic intensity is measured using two pressure-sensing micro- 
phones spaced a distance Ar apart as illustrated in Figure 46. 
The component of the acoustic intensity parallel to the line 
joining the microphones is measured. Let the pressure time 


0 




P = (Pi + P2)/2 
9p/9r = (p 2 - pi)/Ar = Ap/Ar 
SO that Equation (201) for becomes 
= -d/(pAr)) / (p 2 - pi) dt 


(203) 

(204) 

(205) 


p{;ine'pr«“ures"lT“’'^ 


1/ (2pAr)) (pi+p 2 ) /dt (p 2 ~Pi)> 


(206) 


and® d°f/iren®cfi'’'o^^h!®^"^ methods is to construct the 

the pressu« di««ence and ^aJJ^out histories, integrate 

?i5e\^?ouVtTc ^^te/s?ty-.®®®'®^^"* heSCJ?aSrp^Sdic1"US5“?S 

StrSr?Hs.'rEi‘ "•' •■¥•“ ~~ .ssns:? 

terms of the imaaiiarv rfr ^ expressed in 

Gi 2 (f) by, ^«‘agi‘-iary part of cross-^ectral density function 


I (f) = 

r 


Im(Gi 2 (f))/( 27 ifpAr) 


(207) 


Stimate ^She^^cross-spect^ den^i Equation (207) to 
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e.ncti ^ , p?°sT,„:r!'s 

far function H.,f, asl"?alfd wUrra"ch“Mannerrs""* 


J 0 

Hi(f) = IH.(f), e = H^.(0 H^.(n 


( 208 ) 


=P.Pa<^» between the 
i-exatea to the measured Gj 2 (f) by 


‘=P.P2‘" = <^,2(n/IU^(f) H,’'(0 I 


( 209 ) 


cross-spectral ^density ^dS^oted and the measured 
densrty functron is give^t u^n^ier t^h^jse n"e^w%e^s^r?;^^ttL^%^ 




(210) 


the resulting exp^e Js\o^°lea^^ to^”^ taking the square root of 

~ ^^12^^^ {|Hj(f)| iH^Cf)!} (2H) 

the acoustic intensity may be estimated by 

I^(f) = Im{[Gj2(f) ^i2^^^(0]^)/[2nfpArlH^(f)l |H^(f)|] (212) 

averages® o^Tcro"^^^ two ensemble 

product Of the cross-spectral denTV®^ quantity - namely! thi 
addition, the magnitude^ of the !ii!?t^ 

must be measured for each channel. "mentation transfer function 

mentation^^riso\Ti?e^d^ wit^^ ^?h?^n^l instru- 
amplifier, and amplifier in microphone, pre- 

filters and recording devices. If^ont^ucJ? other possible 

and pre-amplifiers which are nhaQA ^ of microphones 

necessary to switch the re™a!ring"lnVt;rent*.ti^ 
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